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PREFACE. 



An abstract of the following investigation was published in Nos. 

(2024-2026) of the Astronomische Nachrichten. But, a few verbal in- 
accuracies, together with numerous typographical errors, being found in 
the work as printed, it was decided to reprint the whole investigation, 
giving the mathematical developments more in detail, and also develop- 
ing the formula for the latitude of the moon, a subject not touched 
upon in that investigation. The author hopes to find leisure for the 
complete development of the perturbations of the moon's motions, by- 
means of the differential equations given in this first chapter. It is 
believed that these equations will give the co-ordinates of the moon in 
a more direct and simple manner than any combination of the general 
differential equations of motion hitherto employed for that purpose. 



Ill 



\ 



Theory of the Moon's Motion. 



INTRODUCTION. 

1. It is now nearly two centuries since the theory of the universal gravitation 
of matter according to the law of the inverse square of the distances was dis- 
covered and subjected to calculation by Newton. The theory when regarded as 
a law of nature may be stated as follows: Every particle of matter in the universe 
attracts every other particle with a force which varies directly as its mass, and 
inversely as the square of the distance between them. 

Mathematicians have attempted to deduce from this law of matter all the 
phenomena attending the motions of the heavenly bodies. For this purpose the 
sun, planets and satellites may be regarded as particles of matter when compared 
with the distances which separate them from the other bodies of the system, the 
masses of the planets being regarded as infinitely small in comparison with that 
of the sun. The principal mathematical consequences resulting from the opera- 
tion of this law of matter may be stated as in the three following theorems : 

I. The orbits of the planets and comets are conic sections in which the sun 
occupies the principal focus; 

II. The radium vector of each planet or comet sweeps over equal area^ in equal 
times; and 

' III. The squares of the times of revolution of the different planets are to each 
other in the same proportion as the cubes of their mean distances from the sun. 

These three laws, which were discovered by Kepler, may be regarded as the 
embodiment of the theory of universal gravitation, and are the foundation of 
physical or mathematical astronomy. They are, however, exact only on the sup- 
position that the masses of the planets are infinitely small in comparison with 
that of the sun, and hence are not applicable to the solar system without some 
modification, because the masses of some of the planets are finite instead of being 
infinitely small. The motions of the planets therefore do not strictly conform to 
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the preceding laws ; for, in obeying their mutual attraction, they must deviate a 
little from the elliptical paths w.hich they would exactly follow if they were 
attracted only by the sun. Not only are the planets disturbed in their motions 
around the sun by reason of their mutual attraction, but the satellites are also 
disturbed in their motions around their primaries by the attraction of the sun and 
by the other planets of the system. The smallness of the planetary masses in 
comparison with that of the sun, however, permits these bodies to conform very 
nearly to the laws of the elliptical motion ; and hence the calculation of their per- 
turbations is not a difficult problem. But in the motions of the satellites about 
their primaries, the principal disturbing body is the sun ; and the eflPect of his 
attraction is too great to be overlooked. Especially is this the case in regard to 
the motion of the moon around the earth, the disturbing force of the sun being 
so great as to cause the transverse axis of her orbit to describe a complete revo- 
lution in about nine years, and also causing the nodes of the orbit to revolve in 
about double that interval. 

2. The mathematical theory of the perturbations of the planets and satellites 
develops the laws according to which the various forces operate in the disturbance 
of each other's movements. They may be classed as follows : 

First The operation of a comparatively large force during a short interval 
of time ; 

Second, The operation of a comparatively small force during a long period 
of time; and 

Third, A very slow change in the mean motions of the diflPerent bodies arising 
from the variations of the elements of their orbits produced by their mutual 
attractions. 

The first class produces the periodic inequalities, which depend solely on the 
mutual distances and configurations of the difierent bodies ; the second class pro- 
duces the inequalities of long period, which may be said to depend on the configu- 
ration of the elements of their motions ; and the third class produces the secular 
inequalities. The second and third classes afiect only the mean motions. 

3. The problem of deducing all the circumstances which afiect the motions of 
the difierent bodies of the system, directly from the principle of universal gravi- 
tation, has, perhaps, more than any other exercised the ingenuity of mathema- 
ticians and astronomers during the last century and a half. The theories of the 
periodic and secular inequalities of the planets were developed in sufficient detail 
for the purposes of astronomy about the close of the last century. Considerable 

progxeaA was also made in the development of the lunar theory by La Place ; but 
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the magnitude of the disturbing forces being so much greater than those which 
disturb the planetary motions, the labor necessary for the construction of a 
perfect lunar theoiy is vastly increased. La Place states that he has determined 
all the inequalities of the first, second and third orders, and the most important 
ones of the fourth order, and has continued the approximation to quantities of 
the fourth order inclusively, and also retained those of the fifth order, which arise 
in the calculation. But much more elaborate theories of the moon's motion were 
published during the first half of the present century, in which the ^proxima- 
tions are carried to terms of a much higher order than waB attempted by 
La Place ; the most complete' and perfect, being that of Plana, in which the 
approximations are carried to terms of the seventh order, was published in the 
year 1832. Pont^coulant also investigated the lunar theory in the fourth 
volume of his Thiorie Analytique du Systtme du Monde, which was published in 
the year 1846, in which he has also carried the approximations to terms of the 
seventh order. The theories of the moon's motion by Plana and Pont:6coulant 
were produced by entirely different methods of development — the one by employ- 
ing the true longitude as the independent variable, and the other using the mean 
longitude for the same purpose — and when reduced to the same form were found 
to be almost identically the same, the coefficients of the inequalities seldom dif- 
fering from each other by so much as one second of arc. It would therefore 
seem that the calculations of both had been correctly made. But perhaps the 
most elaborate and complete investigations on the theory of the moon's motion are 
those of Hansen and Delaunay. Hansen published very elaborate '* Tables de 
la Lune'' in the year 1857; but the analysis by which he obtained his formulae 
has not, so far as I am informed, yet been published. The lunar theory by 
Delaunay was published in two large quarto volumes, which appeared in the 
years 1860 and 1867, and is perhaps the most complete and perfect work on the 
lunar theory that has yet been given to astronomers. 

4. In a valuable and interesting discussion on the apparent inequalities of long 
period in the motion of the moon, published in the American Journal of Science 
and Arts for September, 1870, Professor Newcomb makes the following statement 
in regard to the inequalities of short period in the motion of the moon : 

'* The problem of determining the motion of the moon around the earth under 
the influence of the combined attraction of the sun and planets has, more than 
any other, called forth the efforts of mathematicians and astronomers. Nearly 
every great geometer since Newton has added something to the simplicity or the 
accuracy of the solution, and in our own day we have seen it successfully com- 
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pleted in its simplest form, in which the earth, the moon and the sun are regarded 
as material points, moving under the influence of their mutual attractions. The 
satisfactory solutions are due to the genius of Hansen and of Delaunay. Work- 
ing independently of each other, each using a method of his own invention more 
rigorous than had before been applied, they arrived at expressions for the longi- 
tude of the moon which, being compared, were found to exhibit an average dis- 
crepancy of less than a second of arc. No doubt could remain of the substantial 
correctness of each." 

This statement, coming from so eminent an astronomer, would lead us to 

suppose that the mathematical theory of the moon's motion, in so far as it depends 

on the attraction of the sun and planets, may be regarded as perfect. Four years 

later the same writer, after a thorough and extended examination of the question, 

' assures us of the correctness of his conclusions. 

Hansen has also given in his tables of the m^on two inequalities of long 
period depending on the action of V&nuB, Delaunay, in computing the coef- 
ficients of these two inequalities by his own method, reproduces, very nearly, 
Hansen's coefficient of one of them, but finds the other to be insensible. 

Professor NewcoIib then proceeds to inquire whether we have in either theory 
a satisfactory agreement with observations, and finally arrives at the conclusion 
that the problem of the inequalities of long period in the moon's mean motion is 
really no nearer such a solution as will agree with observation than when it was 
left by La Place; and he is obliged to resort to the operation of irregular and 
extraneous causes, the data for the accurate determination of which are wanting, 
in order to fill up the gap between theory and observation. 

In the Monthly Wotices of the Moyal Astronomical Society for November 14, 
1873, Mr. Airy also enters into a discussion in regard to the inequalities of long 
period in the moon's mean motion, and states his conclusion in these words : " I 
conceive it to be totally impossible that a complete and correct theory can leave 
such large and systematic discordances. And I express my opinion that there is 
some serious defect in the Lunar Theory." 

These statements, coming from such eminent astronomers, are sufficient to 
encourage the belief that the published theories of the inequalities of long period 
in the moon's mean motion possess but very little value. 

6. Let us now inquire whether the inequalities of short period in the moon's 
motion are any more perfectly represented by existing theories than the inequali- 
ties of long period. The theories of Hansen and Delaunay, according to Pro- 
fessor Newcomb, may be regarded as identical in their results ; and they also 
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contain all the important researches of astronomers up to a very recent date. 
The theories of Plana and Pontecoulant may also be regarded as identical in 
their results; the former was published in 1832, and *^ Tables of the Moon,*' 
founded on it, were published by Professor Peirce in 1853. Plana's theory was 
somewhat modified by the theoretical investigations of Hansen ; and some 
empirical corrections, proposed by Airy and Longstreth, were also introduced ; 
BO that the tables in their present shape do not exhibit the results of pure theory 
alone. But if the empirical corrections introduced are legitimate parts of a cor- 
rect lunar theory, and only await the demonstrations of the physical astronomer, 
we may assume that they have a permanent value, and that the tables are perma- 
nently improved by their introduction ; otherwise they possess only a temporary 
value, and will vitiate the tables after a while to the extent of their primitive 
improvement. Now, the ephemeris of the moon, contained in the American 
Sphemeris and Nautical Almanac^ is derived from Peirce's tables of the moon, 
and we have only to consult the volumes of the Naval Observatory at Washington 
in order to learn to what extent the moon deviates from the requirements of 
theory. In this comparison we may suppose that the moon's mean place is cor- 
rectly given by the tables, in which case the observed deviations between theory 
and observation will serve to indicate the magnitude of one or more inequalities 
of short period in the moon's motion which analysis has yet failed to point out. 
We shall therefore find that, in the year 1856, the moon's place differed from the 
computed place between the limits + 10".4 and — 12".4; in 1857 the errors were 
between the limits + 7".7 and — 10".4 '; in 1858 they were between + 14".4: and 
-9".3; in 1861 they were +10".5 and -9".6; in 1862 they were +9'^0 and 
-11".3; in 1871 they were +8".4 and -4".8; and in 1872 they were +14^0 
and — 12'^0. These residuals of course show only the observed deviation of the 
moon from her computed place, and do not necessarily indicate the maximum 
extent to which the theory deviates from the true place of the moon, as cloudy 
weather no doubt frequently prevents observations at times when the deviations 
are at their maximum values. The residuals above given would seem to indicate 
that the mean place of the moon was very closely given by the tables, since the 
positive and negative errors are very nearly equal to each other. According to 
Professor Newcomb, Hansen's tables represented the moon's place at the time of 
their publication and for a few years after much more closely than Peirce's. In 
the year 1862 the errors of Hansen's tables were included within the limits of 
+ 4".0 and - 9".0 ; in 1871 the limits were + 2'^5 and - 16".4 ; while in 1872 
they were — 0".2 and ~16".2; thus indicating a narrower range of periodic 
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variations, and also a much wider range in a progressive mean motion. It should 
be observed, however, that the preceding comparison exhibits only errors of right 
ascension, and it is probable that the errors in longitude would generally be some- 
what larger. 

If we in like manner compare the observed and computed longitudes of the 
moon as given in the Reductions of the Greenwich Lunar Observations from 1831 
to 1851, after correcting for the error of mean longitude, we shall find an oscilla- 
tion between theory and observation amounting to an average value of about 16". 

La Place says, in the Introduction to his Theory of the Moon, that the error 
of the tables formed from his theory will very rarely exceed 32" ; and he also 
states that the astronomer BuRa, by deriving the forms all the arguments from 
theory, and rectifying the coefficients by means of numerous observations, had con- 
structed tables of the moon's motion whose greatest errors were less than 13". 
If this statement of La Place, which was made three quarters of a century ago, 
was borne out by the observations of the moon made at that time, it would 
indicate a degree of perfection in the lunar tables, with respect to the inequali- 
ties of short period, which is scarcely exceeded by the tables in use at the present 
time. 

6. From the preceding comparison it would appear that the lunar theories of 
Hansen and Delaunay, when measured by the criterion to which all physical 
theories must be subjected, are but little, if any, more perfect than the theories 
of Plana and Pontecoulant, which preceded them by twenty years; and we 
also perceive that the lunar theory employed in the Reductions of the Greenwich 
Lunar Observations, and also the tables in use at the close of the last century, in 
so far as the inequalities of short period are concerned, are but little inferior to 
the tables in use at the present time. It is also possible that the superiority 
of recent tables over those in' use half a century ago may arise more from the 
corrections to the mean elements of her motion than from any improvement in 
the theory of her perturbations. It would thus appear either that the theory of 
gravitation had been imperfectly or incorrectly developed, or else that the moon 
undergoes perturbations from the action of other forces than gravitation. 

7. The residuals above given would seem to indicate that our present lunar 
tables, instead of being correct to terms of the seventh order, are really erroneous 
by some of the smaller terms of the third order ; and as the writer had, previous 
to making this comparison, assured himself, by a careful examination of the 
mathematical theory of her motion, that some terms of the third order had been 
overlooked, he does not hesitate to announce the fact to astronomers; and he 
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confidently believes that a correct theory developed to terms of the fifth order 
will be found to represent the motions of the moon, in so far as the inequalities of 
short period are concerned, with far greater precision than any published at the 
present day. 

8. It is true that we have the combined assurance of all the great mathema- 
ticians of the present century who have given especial attention to the subject 
that no inequalities have been overlooked, and that the theory possesses all the 
accuracy claimed for it. To this we would reply that the error is fundamental 
and precedes any development of the perturbing function, simply growing out of 
the latitude of the perigee in the development of the undisturbed elliptical 
motion. If this statement is correct, it is easy to perceive that the fundamental 
error would necessarily ramify the whole development of the perturbing function 
and vitiate more or less all the conclusions deduced from it. Nor do we regard 
the argument that the agreement of so many profound mathematicians, that the 
theory of gravitation when legitimately applied to the moon's motion would pro- 
duce the existing lunar theories, as possessing much weight. Such evidence is 
merely negative in its character, and possesses no importance whatever when 
coming in conflict with positive evidence to the contrary. If any number of per- 
sons develop any mathematical expression by as many diflferent methods, and all 
obtain identically the same result, it is satisfactory proof that the development 
has been correctly made ; but it is no more certain or satisfactory than it would 
be if & singl.e person had made all the developments by as many diflferent methods. 
Besides, if an equation is supposed to possess a certain physical or geometrical 
property, the development of the equation into an infinite series, or in any other 
form whatever, must also possess the same property. But no correct development 
or transformation of an equation can introduce any new element, either physical 
or geometrical. The development must reproduce all the errors or infirmities of 
the function from which it was derived. The assumption that a given mathe- 
matical expression possesses a certain physical or geometrical property should 
therefore be fully justified by a rigorous discussion before such expression is made 
the basis of a physical theory. It is to this point in the lunar theory that we 

« 

would now call attention. 

9. For this purpose let us take the fundamental equations for the latitude 
and reciprocal of the projected radius vector of the moon which are used by 
La Place and Plana as the basis of their respective theories of the moon's 
motion. If we denote the mean distance and eccentricity of the moon's orbit by 
a and e, the longitudes of the perigee and node by w and Q, the inclination of the 
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orbit to the fixed plane by i, and the longitude, latitude and radius vector by 

V, and r, and also put 

y = tan t, « = tan ^, (1) 

« 

we shall have the following equations, in which u denotes the reciprocal of the 

projected radius vector. 

s = y sin(v — S^), (2) 

^=-7T^{|AT7+ecos(i;-6>)}=-i-. (3) 

a(l — r) rcos^ 

La Place and Plana have supposed that equation (3) is the equation of a pro- 
jected ellipse, and we shall now proceed to show that it is not such an equation 
except for the particular case in which the transverse axis lies in the plane of 
projection. 

Equation (2) gives 

l/TT7 = -|/l+y'8in^(v-g3)=l/l+tan^<? = -^. (4) 

cos V 

Whence, cos d = - , (5) 

•|/l + y^sin2(v-a) 

sing= ysin(^~Q) ^ g 

l/l + y«sin^(v-a) 

Now, since u = -, we shall evidently have the maximum and minimum 

rcos^ 

values of u when r and cos d are respectively a minimum and maximum. Now, 

the minimum value of r in an ellipse is r = a(l — e), and equation (5) gives 

cos^= — , when cos^ is a minimum; and the maximum values of r and 

i/T+7" 

COS 6 are a(l+ e) and 1, respectively; we shall therefore find 

maximum value of -w = -^ '— , (7) 

and minimum value of t^ = — — . (8) 

a(l + e) 

If we now substitute the value of \/Y+~?m equation (3), it will become 

•'• u:= ^——Wl+f^in\v-a) + eco^{y-'w)\] (9) 



r cos d a(l — ^) 

and this equation gives 

maximum value oi u = ^ ^ , (10) 

a(l — e^ 
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and minimum value of i^ = — ; . (11) 

a(l + e) ^ ' 

It will thus be seen that while the two minimum values of u agree, the two 

1 
maximum values differ by the quantity -(iey^) very nearly, or by a quantity of 

a 

the third order. We shall soon see that the reason of this coincidence of the 
minimum values arises from the fact that the transverse axis of the orbit is in 
the plane of projection, in which case equation (9) is correct. 

If we now multiply equation (9) by cos Q = , we shall find 

l/TT7^sin> - Q) 

1^ 1 fi+ ecosfc-c.) _ ^^2) 



v — oi) 1 
in^(^r=S^J 



r a{^-&^)\ yT+ysin^ 

In an ellipse the maximum and minimum values of r are a(l + e) and a(l-— 6), 
respectively. If then in equation (12) we suppose that cw = Q, we shall find the 
maximum and minimum values of r to be a(l + e) and a(l— e), as they evidently 
should be ; but if we suppose the ellipse to revolve in its own plane through an 
angle of 90°, oi will become S3 + 90°, and equation (12) will become 

1_ J- fl_^ esin(i;-g3) 1 -^. 

r a{l -^) \ ^l-^fQm\v-Q>) ) ' 
and this equation gives 

maximum value of r = a(l — e?) -«- J 1 * L (14) 

I ^ i/TT?/ 

and minimum value of r = a(l — e^ -^. J 1 + I, 

It will be easily seen that the maxima and minima of r derived from equa- 
tion (13) differ from the true elliptical maxima and minima by quantities of the 
third order. 

Let us now put the first differential coefficient of equation (12) equal to 
nothing, and we shall find, 

dr sin (v — cd) f cos (v — a)) cos (v — Q) sin (v—Q) ^ .^^v 

^^v ~ v'TT7'Z52(^r^~Q) {l+fsm\v-Q)}i "" 

It is evident that this equation cannot be satisfied when v — q) = unless 

v — Q also equal 0° or 90°, which requires that the transverse axis of the orbit 

should be in the line of the nodes or in a line perpendicular to it ; and as this is 

' not the case in nature, it follows that the equation is not applicable to the problem. 



(15) 
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10, Let us now transform the equation of the ellipse in such manner that the 
anomalies, instead of being measured on the plane of the orbit, may be measured 
on the fixed plane which is inclined to the plane of the orbit by the angle i; and 
let the place of the ascending node of the orbit on the fixed plane be denoted by 
Q, while the quantities v, and o) have the same significance as in the last 
article ; also let v' and o)' denote the place of the moon and of the perigee in the 
orbit, and Oq the latitude of the perigee. Eeference to the figure will make the 
geometrical conception clear. 

The equation of the ellipse will evi- 
dently give 

1 1 + 6 cos (V' — fi>') /i fjy. 

r a(l -« er) 

and the required transformation consists in finding the value of cos (v' — o') in 

terms of the angle v — w. Now, if we put v' — o' = (v' —Q,) — (a>' — ^), we shall 

have 

cos (v' — fi>') = cos (v'— ^) cos (o)' — SJ) + sin (v' — Q) sin (a>' — Q). (18) 

We also have 

tan i = y, tan ^ = y sin (v — £J) ; (19) 




cos^ = 



l/l+y^sin^(i;-g3) 



sin^ = 



__ y sin {v — Q) 



cos^o= 



sin 0^ = 



•|/l + >^sin^(fi>-«) 



l/l + )^sin^(v~^) 

y sin (fti — Q) 
l/l+y^sin(a>-a)* 



(20) 



(21) 



Then from the right-angled spherical triangles of the figure it is easy to 

deduce 

cos (v' — SJ) = cos (v — Q) cos ^, sin (v' — Q) = sin ^ cosec i ; (22) 



cos (a>' — S^) = cos (fi> — Q) cos do , sin (a>' — Q) = sin da cosec i. 



(23) 



If we substitute these values in equation (18), it will become 
cos iy' — fi>') = cos (v — Q) cos (oi — Q) cos cos ^o + sin sin ^o cosec* i. (24) 



Substituting the values of sin 0o and cos 0o in this equation, it becomes 

cos {v — Q>) cos (fi> — Q) cos fl + y sin ^ sin (o) — Q) cosec* i 



cos (v' — fi>') = 



Now we have 



and 



l/l+fsin^a) - Q,) 
y sin ^ = y* cos sin (v — Q), 
y* cosec* i = tan*i cosec* i = 1 + tan*t = 1 + y*. 



(25) 

(26) 
(27) 
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If we now substitute the value of ysin^ in equation (25), and then the 
value of -f cosec^ t, we shall obtain, after multiplying by e and putting 

cos C/a = , 



v/l + y2sin2(fi>-«) 

ecos (v' — ctf') == ^ cos^o cos 9 cos (v — S^) cos (ai — S3) 
+ e (1 + j^ cos ^0 cos 9 sin (v — S3) sin (fi> — S2), 

= e (1 + |y^)costfocos^cos(v — o) — Jey*cos^ocos^cos(v +fi>-^2S3). . 



• (28) 



Substituting this value of e cos (v' — fi>') in equation (17), it becomes 

-= {l+e(l+ Jj^) cos^ocos^ cos(v— tt>)— J6y^cos^ocos^cos(v+a>— 2Q)}, (29) 

r a(l— r) 

which is the polar equation of the ellipse, in which the polar angle is measured 
on a plane inclined by an angle whose tangent is y, to the plane of the ellipse 
itself. 

Let us now discuss this equation and see if it meets all the requirements of the 
problem. If we substitute the values of cos^ and cos^o, in the second member 
it will become 

1 __ 1 f 1 e cos (v — co) — \e:f cos (v + cw — 2 S3) 1 /oax 

r a(l-e2)\ i/l+7"sin ^ (^o - S3) i/l + 7' sin^ (v - S3) J 

If we suppose that v = (o, it will give 

a(l— e^ 1 , e + iey^— ie>^cos^(fo — S3) ^ 1 . e(l + y^ sin2(6> — S3) ^ , /oi\ 

r 1 + y^ sm^ (a> — e) 1 +y^siQ*(a>— S3) 

whatever be the relative values of o) and S3. 

The maximum and minimum values of r are also equal to a(l + e) and a(l — e), 
which are their correct values. If we now take the differential coefficient of 
equation (30), we shall find that it may be reduced to the following form : 

a(l - e")— = g(l+y')cosgo8in(t;-^o) .^2) 

r'dv {l + y28in2(t?~S3)}* 

This always becomes equal to nothing when v = co, and when v = co-\- 180°, as 
it manifestly should, since the radius vector is a maximum or minimum at- the 
extremities of the transverse axis. 

If we now divide equation (29) by cos d = - , it will become 

l/l+y'sin2(v-S3) 
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u = -^ = — -— -{l/l + )^siii^(v~g3) + e(l + i}^cos»ocos(v-g2)) 

rcosd a(l—er} I (33) 

— ^ef cos ^ocos (v + tt> — 2S3)}, j 

This is the correct value of u corresponding to an elliptical orbit ; and if we 
suppose that fi> = ^ (in which case cos^o = 1), it will become 

u = -— — — {\/l'\-f8iR^(v—Q,) + e cos (v — o)}, (34) 

a(l— 6^) 

which is the same as equation (3), after substituting the value of |/l+«^. Equa- 
tion (33) gives 

maximum value oi u= -^ ^ , (35) 

a(l — e) 

and minimum value oi u = ; (36) 

a{l + e) 

which are identically the same as before found in equations (7) and (8). 

We thus see that equations (30) and (33) meet all the requirements of the 
problem, and are therefore correct. For the case in which q} = Q>, they immedi- 
ately change to the equations used by La Place and Plana. We also see that 
for the general problem their equation is erroneous by terms of the third order, 
and hence their development of the lunar theory must also be erroneous by terms 
of the third order. 

Having thus shown that the equation which determines the moon's distance 
from the earth, according to the theories of La Place and Plana, is not a conic 
section, except for a particular case, which does not exist in nature, we shall now 
deduce from the general differential equations of the motion of a body acted on 
by the forces of gravitation, the true law of its motion when the positions and 
magnitudes of the different forces are given. 



CHAPTER I. 

GENERAL DIFFERENTIAL EQUATIONS OF THE MOON'S MOTION, WITH THE THEORY OF 
HER ELLIPTICAL MOTION, AND THE VARIATION OF THE ARBITRARY CONSTANTS. 

1. For this purpose let us take the general differential equations of the 
motion of a body which is acted upon by any number of gravitating forces, which 
equations are as follows {M^caniqiLe Cdeate [499], Bowditch's translation) : 

ddx ldQ\ ddy /dQ\ ddz ldQ\ ... 

d^~'\dxy de~'\dy)' 'd^^XdzJ ^^ 

In these equations rr, y and z denote the rectangular co-ordinates, and the 
function Q represents all the forces which act upon the body whose motion is 
required. The time is denoted by t, whose element dt is supposed to be constant. 
We may put equations (A) under a more convenient form for computation in the 
following manner. If we denote the sum of the masses of the moon and earth by 
/£, and all the other forces which act upon the moon by H, the function Q will be 
given by the equation „ 

Q==^^Ii. (1) 

r 

The co-ordinates, x, y and z, of the moon will also be given by the equations 
a: = rcos^cosv, y = rcos^sinv, 2i = rsin^. (2) 

If we substitute these values of x, y, z and Q in equations (A), they will 
become 

ddr — rdv^cos^d — rdd^ u /dli\ .^x > 



2rdrdv cos* <? — 2r* si n ^ cos (? dvdd + r* cos* ddv 



d^ 

2rdrde + r» sin g cos ddv* + r'dde 

de 



ldB\ 

\dv) 



[ddl' 



(4) 



(5) 



(AO 



The integrals of these equations will be the polar co-ordinates r, v and of 
the body whose motion is required. 

If we multiply equation (4) by dt, and integrate, we shall find 

7^ cos' d^=.c-C(^]dt, (6) 

dt J \dv I 

c being an arbitrary constant quantity. 

3 ^ 
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If we now multiply equation (4) by idLn d sinvdi, and equation (5) by 
COS vdt, and take the sum of their products, we shall get, after putting 
2sin'» = l-cos2^ + sin'», 




2rdrdv sin d cos ^ sin v + r^dv^ sin cos d cos v + r^ddv sin ^ cos ^ sin v 
+ r^dvd cos^ sin v + 2rdrdd cos v — i^dvdO sin v + r^ddO cos v 
— T^dvdd sin^5 sin v 

= — d^J ( — jtan^sinv + ( — jcosv L 
\\dv/ \d0J J 



Equation (7) gives by integration, 

1 

— {T^dv sin d cos ^ sin v + r^ cZ^ cos v} 

dt 

= c'— \dt\ I — Itanflsinv + ( — |co8v i 



. (7) 



(8) 



In like manner, if we multiply equation (4) by — tan d cos vdt, and (5) by 
sin vdt, we shall get, by integrating the sum of their products, 

1 



dt 



{r^dO sin v — T^dv sin cos ^ cos v} 

= c"+ Cdt\l — Itan^cosv— (-— -|sinv I; | 

J \\dvj \do) y J 



(9) 



c' and c" being arbitrary constant quantities. 

If we now multiply equations (A') by 7^{dv cos ^ sin v + d0 sin ^ cos v} , 

dr +2rc?vcos^cosi;, and c?r sin ^ cos v + 2d^ cos ^ cos v, respectively, and 

cos fl 

take the sum of the products, we shall obtain 

Sr^drdO^ cos ^ cos v + 27^d0dd0 cos ^ cos v — 7^d0^ sin cos v 
— 7^d0Hv cos ^ sin V + Zr^drd'i^ cos* tf cos v + 2ii^dvddv cos*^ cos v 
—Sr^dv^dO sin cos* ^ cos v — r^drP cos' ^ sin v + 2rd7^dv cos^ sin v 
+ r^ddrdv cos 0Bmv + r^drddv cos ^ sin v — 7^drdvd0 sin ^ sin v 
+ r^drdv^ cos ^ cos v + 2rdr^d0 sin ^ cos v + r^ddrd sin tf cos v 
+ 7^drdd0 sin ^ cos v + r^drd0^ cos ^ cos v — r^drdvd sin ^ sin v 
+ /i {cos ^ sin vdv + sin ^ cos vd0} 

/dli\ 



dfi 



— — r*{dvcos^sinv+rf^sin^cosv}( — ]-'{dr-^^+2rdvcoQ0co&v}[- 

\dr I cos^ \dv ) 

— {dr sin ^ cos v + 2rd0 cos ^ cos v] (-— ■ ) 
Multiplying equations (AQ respectively by r*{<i^sin^sinv — dvcos^coBv}, 



(10) 
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cos V 
— dr + 2rdv cos Q sin i?, and dr sin ^ sin v + 2rd9 cos 9 sin v, the sum of their 

cos 5 
products will give 

Sr^drdd^ cos fl sin v + 27^d0ddO cos fl sin v — r'dfl^ sin ^ sin v ^ 

+ T^dvdO^ cos 5 cos V + Sr^drdv^ cos*/? sin v 4- 27^dvddv cos¥ sin t? 
1 — Sr^dv'^dd Qmd cos^fl sin v + r^d-y* cos*^ cos v — 2rc?r^cZv cos ^ cos v 
^ I — r^ddrdv cos cos v + rdrdvdO sin ^ cos v + r^drdv^ cos ^ sin v 
— T^drddv cos ^ cos t? + 2rdr^dd sin /? sin v + T^ddrd d sin tf sin i? 
4- r^drddO sin /? sin v + r^drdff^ cos fl sin t? + r^drdvdO sin ^ cos v ► • (1 1) 

+ ;i{(i/? sin ^ sin v — cZv cos 5 cos v} 

^r^fcZv cos ^cos V — dfl sin sin v} ( — ) + {dr — 2rdvcos ^sin v} ( — ) 

\dr / cos d \ c^'y / 

— (drsin^ sinv + 2r(i^cos^sinv}( — ) 

If we now multiply equations (A') respectively by —T^ddcosO, 2Tdvsixi0, 
and — dr cos d 4- 2rdO sin d, the sum of their products will give 



if 

de ] 



V 



—2rdr^dd cos ^ — r^ddrdd cos tf — r^drddd cos ^ + r^drdd^ sin ^ 
+ Sr^drd^^sin » + r^dv^^^ ^.^gS^ _|_ ^r^ddddd sin » + r^dO^ cos» 
+ Sr^drdt?^ sin^ cos^^ + 27^dvddv sin ^ cos^fl — 2r^dv^dd sin^O cos 
— )M cos ^ dd 

^dd cos <? /— ] - 2rdv sin <? (— \ + {dr cos 0- 2rde sin »} /— \ 



(12) 



Lastly, if we multiply equations (A') respectively by 2dr, 2dv and 2dd, tie 
sum of their products will give 



1 J 2drddr + ^-drdv^ cosW + 2rdrdff' + 27'dedde 1 , « ^ "^ 
de [+ 2r'dvddv cos» 6 - 2r'dv^de sin ^ cos ^ J '^ r» 



= -2 



'dM\ 
dr) 



dr + 



(f)--(f)-} 



(13) 



If we now take the integrals of equations (10), (11), (12) and (13), we shall obtain 

— {r*cZ(?^ cos^ cos V +r'cZv^ cos^tf cost? +r^cZr dvcos^ sin V +7^(ir<i^ sin^cosv} 
dr 



—ficoQdcoav—f 
= — r J 7^{(ivcos^sinv+<i^sin^cos 



v}l — )+{c?r k2rc?vcos^cosv 

\dr/ 



cos^ 



+ {dr sin ^ cos v + 2rdd cos cost? 



>(f) 

m 1 
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—- {r*ci^co8(?sini; + r'di^ cos? Oa.nv — i'drdv cos dco&v +r'drdd Bin Oainv} 
or 



—ft cos Oainv —f 
— r J r*{d(? sin^sinw —dv coB<? cosw} [-r— j+{2rrft) coaO sinv a'^^'H "T" I 

+ {draiad aiav + 2rd6 cosd smv}l^-\ \ 

\d0J) J J 



L(i5) 



d«« 



{r'djO* sin <? + r»(i'»« sin ^ cob* ^ - I'drdS cos ^} - ;i sin <? -/" 



.f jr» (i<? cos <? /^ - 2rdt> sin <? /^ -{2r<i<? sin tf - dr cos <?} /^ I 



; (16) 






a 



= -2 



/{(f) 



c?r+ 




dv + 




(17) 



/i /'> /'' ^'^d -- being axbitrary constant quantities to complete the integrals. 

a 

Equations (6), (8), (9), (14), (15), (16) and (17) are the polar equivalents of 

equations (P), M^cxmique CHeste [572]. 

If we now multiply equation (8) by cos v, and equation (9) by sin v, the sum 
of their products will give 



dd __ c^cosi; + c^^sinv 



r^ 



+ 



^/*{(f) 



tan 5 cos V— 




sinv L 



> • 



— -J— C<ii< i-r^^tan^sinv+ 




/dB\ 

[del 



cosv 



1 



(18) 



Equation (6) gives 



dv 



'd 



dt T^ cos^ 7^ cos^ 6J\dv 



\dt. 



(19) 



If we now multiply equations (14) and (15) by sinv and ~cosv respectively, 

dv 
and divide the sum of their products by r*cos^— -, we shall get 

dt 
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dr f sin v — /' cos v dt^ 
dt r^cos^ ' ' dv 



cos vdt 
r^cmOdv 



f 



Bmvdt 
T^cosOdv 



f 



7^{ddBinOBmv — dv cos ^ cos v} | — | 

\drj 

+ \ 2rdv cos ^ sin v — dr — - I ( — ) 
( cos^J \dv/ 

+ {2rd0 cos ^ sin v + dr sin sin v} I — ) 



7^{d0 sin ^ cos V + dv cos sin v} 



\ 2rdv cos ^ cos V + dr 1 (-^ 

( cos^ J \dv 



> . <20) 




+ {2rd ^ cos 5 cos t? +d!r sin ^ cos 



"'(^J 



The integrals of equations (18), (19) and (20) will give the latitude, longitude 
and radius vector of the moon at any time t To determine the integrals of these 
equations it is however necessary to know the values of the constant quantities 
^> ^\ <^"> fj /'> /" ^^^ ^ which were introduced by the integrations, when the 
function JR is equal to nothing, and we shall now proceed to determine them. 

2. For this purpose we shall suppose that iJ=0, which gives (-i— ) = 0> 
(^\=0, and /— \= 0; and equations (6), (8), (9), (14), (15), (16) and (17) 



will become 



r^cos^tf-r- = c: 
dt 

9 ' n n ' dv m dO , 
IT sm cos u sm V \- ir cos v— - = <r : 

dt dt 

9 ' d0 o . a a dv M, 

r sm V r^ sm ^ cos cos v — = c'\ 

dt dt 



(21) 
(22) 



(23) 



1 f r'c?^ cos ^ cos v + r'civ^ cos* ^ cos V 1 ^ >. /o^x 

— \ ^ — /££ cos ^ cos V =/; (24) 

d^ { + r^drdv cos tf sinv + r^drd0 sin ^ cost? J 

1 fr'dtf^cos^sinv + r'dv^cos^^sint; 1 a - xi /ok\ 

— \ [-'/jLQoa0amv=f'; (20) 

J^ I — r^drdv cos ^ cos v + r^drd0 sin ^ sin v J 



(ft^ 



d^ 



{r'dO^ sin + r^dv" sin cos^^ - i^drdJd Qoa0}^/Jt sin =f' ; 



(26) 
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— {dr" ■Vr'dd^ +r'dv' QO^H) --"^ + f^ =. 0. (27) 

df T a^ 



Equations (18), (19) and (20) also become 

dO __ c^cosv + c^^siny ^ 
di^ ? ' 

dv c 



dt r'coQ^d' 

dr _ f^mv—f'co^v dt^ 
dt T^coQd ' dv 



(28) 



(29) 



(30) 



If we substitute the value of dt given by (29) in equation (28) it will become 

= — cos vdv + — sin vdv . (31) 



cos^tf c c 

Equation (31) gives by integration 

tan d = — sin v — — cos v. (32) 



c' . c" 



When the moon is at the node v = Q, and ^ = 0, therefore equation (32) 
will give 

= — sin SJ — — cos SJ; (33) 

c c 

therefore tan Q = — . (34) 

c' 

rJfl 

When 5 is a maximum it is equal to the inclination of the orbit ; and — = ; 

dv ^^ 
e(juations (22) and (23) will therefore give, by substituting the value — , 

dt 
c' = c tan sin v, c" = — c tan cos v. (36) 

If we square these equations the sum of their squares will give 

tan'g= ^ ^^ =f. (36) 

The quantities c, c' and c" therefore determine the place of the node and the 
inclination of the orbit ; and if we now put 

— = y cos Q, — = y sin Q,, (37) 

C 
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equation (32) will become 

tan^ = y{sint? cos Q — cost? sin Si} = T'sin {v — Q), (38) 

which is the same as equation (2), Int., and gives • 



8in&== ^^ ^ , cos5 = 



Vl+fsm^{v-Si) i/l+y*8in*(v-gi) 



(39) 



If we now square the values of c, d and c", and put their sum equal to A*, we 
shall find 

A3 = ^ + g/2 + ^m = -^{r^ (^^Qd'i? + T'dS'). (40) 



This equation gives 

di? " r»* 

If we substitute this in equation (27) it will become 



(41) 



At the extremities of' the transverse axis, dr = 0, and equation (42) will 

give 

[17* — 2fjar = — ah. (43) 

I ^ 

This equation gives r = a + a\/l-~ — , (44) 

^ /la 

I p 

and r = a — a\/l--— . (45) 

The sum of these values of r is equal to the transverse axis of the orbit, and 

their difference is the double of the eccentricity. Therefore 2a denotes the trans- 

1 p 

verse axis of the orbit, and \/l denotes the ratio of the eccentricity to the 

^ fia 

semitransverse axis or mean distance. If we put 

e = Jl^; (46) 

equations (44) and (45) will give 

r = a(l±e), (47) 

which are the values of r at the extremities of the transverse axis. 
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If we now denote the longitude of the moon, when r is a minimum, by o), and 
also put dr = 0, and v = (o, in equations (24) and (25), they will become 

\—{7^de^coQd +r^dv^ coQ^d} - ;icos^ I co8<m =/; (48) 

f JL{r8d^ cos + T^dv" cos'^} - fi cos ol sin oi =/'. (49)^ 

These two equations give 

tanfi> = 4- (60) 

/ 

If we now take the sum of the squares of the values of /, /' and /", given 
by equations (24-26), we shall obtain 



-^{r'idd^ + dv'Qo^df + r'dr^dO' + dv^ cos»<?)} 



or J 



=/*+/'»+/"* -A*. (61) 



de 
But equation (41) gives 



d0* + dv'coB'd = h^. (62) 

Substituting this in equation (51) it becomes 

^ + h^ - ^ = f +/'* +/"* - fi?. (63) 

r* d^ T 

If in this equation we put dr = 0, r = a (1 dt e), and h^ = ajw (1 — eF), it will 
become 

f +r +/"' = i"V. (64) 

*! cos udt) 
If we now substitute dt= , r = a(l — c), i;==fi>, <? = ^o and<ir = 0, 

. c 

in equations (24), (26) and (26), they will become 

f-71 A 7^-71 + — ^]- /A COS do COB 0) =/, (65) 

La (1 — e) ( cos* ffodtr cos^^oJ J 
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But equation (38) gives, when v = a}, 

dd = y cos^ 00 cos (fi> — Q) dv, (58) 

dd^ 
-whence we get --— — — - = y* cos^(fi> — ft), (69) 

av^ 008*^0 
and — ^ = 1 + y> sin^C^; - Si), (60) 

COS^ do 

whence we get ^J^ + _1_ = 1 + y., (61) 

and equations (55-57) become 

l^^^^-/i|cos<?ocosa> =/, • (62) 

ia(l— 6) J 

I ^7T^ - /i I cos <?o sin a> =/', (63) 

(a(l— 6) J 

But equations (36), (40) and (46) give 

(?{l+f)=(? + c'^ + c''^ = h^^a/i{l-^), " (65) 

therefore we get — ^ = /i{l + e). (66) 

a{l—e) 



Substituting this value in equations (62-64), we get finally, 

/ = /ie cos 00 cos o), 

/' = /ie cos ^0 sin (o, y (67) 

/" = fjte sin ^0- 



1 



Therefore the quantities /, /' and /" denote the product of the sum of the 
masses of the moon and earth into the co-ordinates of the centre of the orbit 
when referred t6 the focus as the origin. 

3. Having thus found the values of the constant quantities introduced by the 
integrations, if we now substitute them in the difierential equations of the co- 
ordinates- r, V and 0, we shall obtain, by means of another integration, the values 
of these co-ordinates in terms of the time. But as we have already found the 
value of in terms of v, in equation (38), we shall also find the values of r and t 
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in terms of v. We shall then, by inverting the formulas, be able to find the 
values of r, v and in terms of the time t. 

Equation (66) gives 



^^Vf^ajl-^^ . (68) 

VT+7 

If we substitute this value in equation (29) it will become 

^=-Vl±Z_.r«cos»«. (69) 

And if we also substitute these values, and also the values of / and /', in 
equation (30), it will give 

dr 6(l+}^)coso s/j . / X 7 /t7n\ 

_ — _3 Li . cos'* tf sm (v — o}) dv. (70) 

r a(l— 6^ 
But we have . cos'^ = — ■ : (71) 



which, being substituted in equation (70), gives 

dr __ e(l +/)co8^osin(i;— a>)c?i; 
7^a(l-e^{l + y2sin2(v-Q)}** 



(72) 



This is identically the same as equation (39), Int. ; and hence equation (72) 
will give, by integration, 



r a(l-e^l 



^ g(l +iy*)co8<?ocos(v — o) — jef COS 00 COS {v + at— 2Si) ) . /i,o\ 

l/l + fain'{v-Q,) )' 



which is the same as equation (36), Int., obtained from purely geometrical con- 
siderations, and is therefore the equation of an ellipse. 

If we now substitute this value of r in equation (69) it will become 

dt^ at(l-e^^V^T+7co8» g ,^^. 

dv /- f 1 e(l + ^y') cos <?o cos (i> — <o) — ^e}^co8<?oC08(p + a»— 2Q) | 
^^l ^l+fain'iv-^ J, 

X 

This equation expresses the rigorous relation which exists between the element 
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of time dt and the differential of the moon's motion in longitude corresponding 
to any part of the orbit ; and its integral will give the true time ^, which is 
required for the moon to pass through any arc of longitude which is denoted by v. 
It is not, however, readily integrated in its present form; but we may develop the 
second member into an infinite series arranged according to the ascending powers 
of e and y, which, in the lunar theory, are numerically small quantities, and carry 
the approximations to any degree of accuracy which may be necessary. The 
different terms of the equation when thus developed can then be integrated 
separately, without any analytical difficulty. We shall now attend to this trans- 
formation of equation (74), and shall carry the approximation to terms of the 
seventh order of magnitude depending on the eccentricity and inclination of 
the orbit. This degree of approximation is greater than is necessary in the 
theory of the moon's motion ; but as the same formula may be applied to the 
motions of the planets, it was thought best to give it all needful extension for 
that purpose. 

4. If we develop the variable part of the second member of equation (74) by 
the binomial theorem, it will become as follows : 



jnf^ , e (1 + \-f) cos Oq cos (v — w)— \ef cos Op cos (v + q> — 2Q) 1 ""* 






yi+fsin^iv-a) 



ta n e (1 + iy^ coa d„ COS (v — a)) — iey* COS 6tCoa(v + w — 2Q) 

= cos P — J . — ^ ^ ^ ; - 

{I + f Bin' (v - Q,)}i 



, „ {e (1 + ^y*) cos 00 cos (v — w) — ^ey^cos <?o cos (t> + a> — 2Q)Y 
' {I + f Bin^v - Si)y 

_ J, {e(l +^y')co8tftCOs(p— a)) — |ey*co8<?ocos(t) + a> — 2Q)}* 

{l+fBm'iv-gi)}i 

,-{«(! + if) cos 0t cos (v — to) — \ef cos 6^ cos [v + w — 2S^) }* 

{l+y*sm*(v-Q)}» 

_ /. {e (1 + \f) cos ^0 cos (v — a> ) — ^ey^cos^ocos {v + a> — 2^)}' 

{l-^fBin'{v-Q)}i 

, ,_ {e (1 + \-f) cos (?o cos {v — (o ) — \efcoBad cos (y + a> — 2Q)}* 

{l + /sin*(u-gj)}* 

* _ jj {e (1 + \-f) cos di, c os (t) — a>) — ^ef cospg cos (« + a» — 2Q) y 

{1 + f sin* (v - g^)} t 



(75) 
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We shall now develop each term of this equation separately. Since 

coaO = , we shall obtain 

l/H-y*8in»<»-gJ) 

C0S»^={H-y«sin»(»-a)}-* = l-y*sin»(v-g3) + /sin^(w-gJ)-/sin«(v-gJ) 



+ i/ (1 - f/) cos4 (w - ft) + ^/cos 6 (t; - Q) 



} (76) 



We also have 



cos ^, = {1 + y» sin* (4u - ft)}- J 

=i-i/+Ay*-is%y'+iy*(i-}y*+^/)co82(a,-ft) 

+ A/ (1 -if) C034 (a> - ft) + yfr/ cos 6 (tt» - ft) 



}• (77) 



We shall therefore obtain the following values of the different terms of the 
second member of equation (75) : 

eQ- + ^7*) OPS 00 co6{v — to) — ^ey* cos^o cos (y + a> — 2ft) 

{l+y«8in*(v-ft)}t 

+ |cy*(l- Z+My*) cos(3v - o» - 2ft) + ley* (1 - f+Hy*) cos {v -3«>+2ft) 

- xf^ (1 - f 7*) cos (3v + tt» - 4ft) + tJ7«/(1 - f y*) cos (v + 3<w - 4ft) 

+ TyV«y*(l-f?^co8(6u-a»-4ft)+Tf^(l-f>^cos(t;-5a>+4ft) i J (78) 

~Tf5Vcos(6v+a>— 6ft)+yfy«/co8(v+5a»— 6ft)+T||jej^cos(7v— <w— 6ft) 

+ yMr^y* cos (6t; - 3tt» - 2ft) + TT^^cos (3w - 6<w + 2ft) 

+ Ti^fT^y* cos (v — 7«> + 6ft) 



{6(1 + ^7^ cos (i>-a))-^cos(» + «-2ft)}*cos»^,{H-y*sin*(t>-ft)}-* 

+ ieV (1 - y*) cos (2w - 4a» + 2ft) + T5>2«V* cos (2?; - 6a> + 4ft) 
+ i^^*cos 4 (v - a») + ^eV* cos (6t) - 2<» - 4ft) 
+ ieV (1 - y*) cos (4» - 2« - 2ft) 



i ; (79) 
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{e(l +irO C08(w-a») -^cos(v + a»-2S2)}'cos»»o{l +y*8in='(t>-^)}-4 

+ ih^^ cos5(v - (o) + AeV(l -7^ cos {Zv-w- 2Q) 

-^f^* cos(3v + fi> -4Q) + ^^^ (1 -/O cos (v - 3a> + 2Q) 

+ ^ey (1 - f) cos (Sv - 3a> - 2g2) + -^r '(1 - y") cos (3v - 6ai + 2g2) 

- -yf^eV* cos (3v + fi> - 4Q) + iri^eV cos (z; + 3fi> - 4Q) 

+ yf^eV^cos (v - 5a> + 4Q) + V^e^* cos {Zv - 7fl> + 4gJ) 

+ ^eV*cos (7n^ - 3a> - 4Q) 



; (80) 



{e (1 + Jy^ cos (v - a>) - -^cos (v + oi - 2Q)}*cos*5o{l + y^sin (v - Q)}-« 

= |6*(l--i>^) + i6*(l--J}^cos2(v-ft>) + ie*(l---^7^cos4(v-iw) 

+ A«Vcos2(i;~g2)+^eVcos(2v-4ft) + 2Q)+^Vcos(4v-2fi>-2Q) j^ ; (81) 
+ ^eVcos (6v ~4fi> - 2Q) + ^^^eVcos (4v - 6fi> + 2Q) 






{6(1 +i>^)cos (v -fi>) -i6}^cos(v + o> -2Q)}'^cos5o*{l + y*sin2(v -i^)}-* 



= fe* (1 — ^y^ cos (t; — ctf ) + ^^e* (1 — ^f) cos 3 (1; — cw ) 
+ 3^e* (1 — Jy^ C085 (v — cw) + ^e'y^cos (3v - 
+ yI? e^y^cos (v — 3a> + 2Q) + i|^e*y^cos (5v ■ 
+ ^«Vcos (3v - 5a> + 2^3) + i-5t«Vco8 (7n^ 
+ xbe^^cos (5v - 7<w + 2Q) 



ai-2Q) 
3ft) -2Q) 
•6ft>-2S3) 






(82) 



{e(l + \-f)co& {v — ai) — ley' cos (v + a» — 2g3)}*coa*<?o 

{l+fsin^v-Q)}* 

= ^y* + ^e*cos2(t> — o)) + -]^e'cos4(v — «>) +-5^cos6(t> — a») 



I 



(83) 



{6(1 + -^y^ cos {v — ft)) —^ey^ cos (?; + ft> — 2Q)}^cos^^o 

= "ll-e^ cos {v — a)) + ^^e^ cos3(v— ft>)+-^e^cos5 (v — ft>) + -g^e^cos 7('y 



-S3) J 



. (84) 



If we now multiply equation (78) by -2, (79) by +3, (80) by -4, (81) by 
+ 5, (82) by — 6, (83) by + 7, and (84) by — 8 ; and then add the products to 
equation (76), we shall get the following development of equation (74) : 
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Vl* 



dt 



a*(l-e*)*l/T+7 dv 

= 1 - iy» + |e^ + 1/ + J,/^ - A/ + f|e« - feV + ^J^ - Hey 

-{e'-^+^ej'H J^+ifeV-T^'+^^-Hey} co8 3(r -0.) 
+ {|«* -• A«V + A«V* + H«*} cos 4 (t; - «) + ^ cos 6 (v - tt») 

- {|c» + -Je^ — ^ey + ^eV*} cos 5 (w - «») - -Je^ cos 7 (t> — <u) 

+ ^/{l - fy*} cos 4 (i> - Q) + ^y« cos 6 (v - g3) 

- dey* — |e/ + }ey - feV + Hk^ + tfey } cos (3v - a> - 2Q) 

- {Jey* - iey* + i^y - iey + ^^^Zff«}^ + if ey } cos (v - 3a> + 2^5) 
+ { A^y' - xfiE^y* + tIW'*} COB (3v + «» - 4^3) 

- { A^y' - TfT«y* + TJW''} cos (v + 3 «> - 4a) 

_ {||ey< - ^ey« + ^e»j,<} cos (5u - a» - 4Q) 

- {ifW - rfsV + xfre^'} cos (v - 5a> + 4g3) 

+ ^f^ey* cos (5t> + <a — 6Q) — -^^-^ cos (v + 5<w — 6Q) 

— ^ey» cos (7n< — o) — 6Q) — ■j|35€)^ cos (v — 7o> + 6Q) 

+ {fey - fey + ifey} cos (2» - 4a> +2Q) 

+ {fey - fey + H «y } cos (4v -2a>- 2Q) 

+ ^je^< cos (2w - 6a> + 4Q) + ^r€y cos (6u - 2fi» - 4^2) 

- {fey -• ley + Trh^» + VV«V} cos (3iJ - 5«, + 2gJ) 

- {fey -■ fey + ?SW + H«V} cos (5 V - 3«» - 2Q) 

+ M«V cos (61; - 4tt» - 2Q) + ^y cos (4n< - 6«> + 2^2) 

- ^^«y cos (3n« - 7<w + 4Q) - T^^* cos (7u - 3«> - 4Q) 

- Hey cos (7n< - 6a» - 2SS) - -ifey cos (5u - 7<w + 2^3) 



(85) 
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Now we have 

(l-c^ivTT?=l-|€« + iy» + |e*-i/-|«V+^V+AcV+A«*+TV/- (86) 

Multiplying equation (85) by this value of (1 — e^)ii/TT^, it becomes 

+ f€^{l + i«' + i^-i/}cos2(v-iw) 

+ {f^* + A^* + 1^} cos 4 (v — tt>) + -^e^cos 6 (v — co) 
"{1^ + A^^ + M^} cos5(v — o>) -|e^cos7(v — tt>) 

.+iy'{i-iy'-}^+Ay*-i«*+l^}cos2(t;-gj) 

+ iy*{l-r'~f^}cos4(v-a)+^y*cos6(v-a) 

-f^{l-iy'-i^+Hy'-iV' + i^}cos(3v-fi>-2gS) 

+ A^{l~)^~i^}co8(3v+^~4Q)-^6y*{l-)^-i6'}cos(v+3fi>---4Q^ 

- ifey^ll—y' — 6*}cos(5v — tt> — 4Q) + yfg^*co8(5v + o> — 6Q) 
— -^^y* {1 — y* — e*} cos(v — 6fi> + 4Q) — -yi^Vcos (v + 6fi> — 6gJ) 
+ |eV {1 -^f-ie"} co8(2v ~4fi> + 2Q) - /lA^ey^cos (7t> - fi> - 6Q) 
+ fcV{l - ^r^-i^} cos(4v - 2fi> - 2Q) -^f^«co8(t; - 7fi> + 6gJ) 

- {|«V'- A^* + t!t^} co8(3i;-5fi> + 2Q) +^6Vco8(2v-6a> +4g3) 

- {|«V - A«V + :bS^^*} cos (5v ~ 3fi> - 2Q) + ^e^* cos {6v -2(0- 4:Q) 

- ^^^* cos (3v - 7ai + 4Q) - ^^* cos (7v - 3a> - 4g3) 
+ if^V cos (6v - 4fi> - 2Q) + ^V cos (4v - 6fi> + 2gS) 

- lieV cos (7v - 6ai - 2g3) - ifeV cos (6v - 7fi> + 2^3) 



If we now put n = 5-^, (88) 

a* 

equation (87) will give, by integration, 



.(87) 
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nt = v-2e(l— ^y* + jij^ + ■fyy') sin {v -co) 
+ |e» (1 + i€* + ^ - i/) sin 2 (v - a;) 

- (ie» + ie» + ^JW* + A«' - iAr«V - TSW) sin 3 (v - «») 
+ {^e* + ^ + Tft«V} sin 4 (iJ - «.) + Ti7«' sin 6 (v - fl>) 
-{^ + ^ + -^*}sin5{v-a))—^'Bm7{v-a>) 

- iy'ii - iy* - 1** + Ay* - K + l«y} sin 2 (^ - «) 

+ T^/ {1 - y* - F} sin 4 (« - SJ) + Ti^r/ sin 6 (t> - ft) 
-AyMl-y*-F}8in4(«»-ft)-Ti3r/sin6(o>-ft) 

- i«y*a - iy* - i^!* + H/ - T^ + i«V} sii^ (3« - «> - 2ft) 
-i«y'{i-iy'-i«*+Hy*-A«*+i<V}sin(«-3a»+2ft) ).. (89) 

+ ^ey^fl-y'-ie^lsinCSt; + a» -4ft) — 8lj«j'*{l-y*-^}sin(t;+3«)-4ft) 
— •^«}'*{ 1 — y* — e*} sin (5» — a» — 4ft) + ^IfC)^ sin (5v + o) — 6ft) 
— ^€y*{l — y* — e*} sin (■» — 5a) + 4ft) — •j^^cy'sin (v + 5tt» — 6ft) 
+ A«V { 1 - iy* - i«*} sin (2i; - 4a> + 2ft) - :p|7ey« sin (7v - a» - 6ft) 
+ A«V'{l-iy*-l«'}sin(4v-2fi»-2ft)-yf^8in(i>-7tt» + 6ft) 

- {ieV - tV«V + A«y*} sin (3» - 5«) + 2ft) + -^^y sin (2» - 6a> + 4ft) 

- {i«V -T^' + Trf^Vlsin (6i; -3o> - 2ft) + V^«y sin (6i;-2a> -4ft) 

- jfyeV* sin (3u - 7a» +4ft) - yfgeV* sin (7t) - 3<o - 4ft) 
+ ^V sin (6w - 4a> - 2ft) + ^^eV sin (4i; - 6«» + 2ft) 

- ^^* sin (7v - 5«> - 2ft) - ^eV sin (6» - 7<w + 2ft) 



In this equation nt denotes the mean longitude of the moon, and v the true 
longitude ; and the constant introduced by the integrations has been made to 
satisfy the condition that the mean and triie longitudes shall be equal to each 
other at the extremities of the transverse axis. By inverting this series we may 
obtain the value of v in terms of nL 

5. For this purpose we shall observe that»if we change the sign of all the 
terms of the second member except the first, and at the same time change v into 
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nt in those terms, and calling the sum of the terms thus changed /(nj), we shall 
have, according to the theorem of La Gkanok, 



V = 






ndt 



n^d^ 



5T 



cPfinty ] 



rv'd^ 



. .dyinty d'>f{nty 

^ 12 ,4 7^ + T1T5 --^^ + 



nW 



TlW 






(90) 



We shall therefore have 



/(7i<) = 2e {1 - T^/ + Tfir^V* + A/} sin (n< - fl.) 

' -|«'{l + ie' + -Ar«*-iy'}8in2(nt-a.) 
+ (i** + F + TtW* + T^' - A«y - IJW) 8i° 3 (»« - a,) 
~ {^«* + ^r«* + A*V*} sin 4 (n< — o) — yJj^ sin 6 (n« — o)) 
+ { A^ + tW + A^*} sin 5 (n« — a>) + ^c' sin^ (w< — tt») 

- iy* {1 - i/ - F + Ay' - i** + l«V} sin 2 (n« - Q) 
+ iy* { 1 - iy* - F + A/ - i«* + l«V} sin 2 («» - Q) 

— ^y* { 1 - y» - 1**} sin 4 (n< - Q) — i^/ sin 6 (n« - ^) 
+ A/ { 1 - y* - F} sin 4 (a» - a) + yi^y* sin 6 («. - ft) 

+ i«y* { 1 - iy* - F + Hy* - A«* + i«V} sin (3n« - a* - 2Q) 

+ ^^^{1 -iy*-i«' + Hy*-A«' + i«V} sin(ne-3«* + 2fl) 
— 5^*{l-y»-i6*}sin(3n«+tt>-4Q)+^'{l-y*-ie'}sin(n«+3«>-4a) 
+ A«y' {^—f-^] sin (6n«-o) -4^5) --jf^ey'sin (5n« + tt>-6g3) 
+ A«y' {1 - y* - «*) sin (n< - 5«> + 4^2) + Trhey« sin (n« + 5«> - 6Q) 

- A^V' { 1 - iy* - F} sin (2n« - 4«> + 2^3) + jfi^ey" sin (7n« - o - 6g^) 

- A^V' {1 - iy* -• K} sin (47i< - 2a> - 2gJ) + ^ey« sin (n« - 7<w + 6^2) 
+ {i«V - tW* + Tfir«y*} sin(3ni;-5«> + 2S3) - -^eV* sin(2n<-6<w+4S^) 

+ {i«V - Ae»y* + sh^^] sin (5n« - 3«> + 2a) - ^e^* sin (6n« - 2a»-4gJ) 
+ T^e^* sin (3n« - 7ft» + 4Q) + yff e»y* sin (Int — Za>- 4^) 
— A^V sin (6n« - 4a> - 2Q) - ^V sin (4n« - 6<i> + 2^) 
+ A«V sin (77i< - 6tt» - 2£J) + ^eV sin {6nt - 7a» + 2^) 



• (91) 
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From this equation we get the following values of the different powers of 
f{nC)f which are all correct to terms of the seventh order : 



- {f e» +^6^ + iey* + ^je^ - i|€y - ie/} cos (»< - a>) 

+ {!«'- tV* + iey^ - A«' -I'i^VeV - i/} cos 3 (n< - «.) 

- {fie' + T^ + A«V'} cos 4 (n< - «») — H||e» cos 6 (n« - 0.) 
+ {A^" + ¥ + A^} cos 5 {nt -<a) + f^e'cos 7 («.<*- ai) 

+ {H«V-H«V'-if«V + TiTf/} {cos2(n«-Q) + cos2(«>-.a)} 

- {hy* - -sV/ - A«V'} {cos 4 {nt - g3) + cos 4 (a. - g^)} 

— Yj^y^{cos6(n^ — Q) +cos6(fi> — Q)} — yf^^y* cos (Sti^ + 3tti — 6^) 

+ {i^-i^'-f^ + AeV*-Tl^eV + ^ey«}cos(n^-3a> + 2Q) 

+ {i^* - i9^ -- i^^^] {cos (n^ - 5fi> + 4Q) + cos {bnt - a> - 4Q) } 

-{tt^'--tt^~i^V+Tky'}{co8(2ri<-4a>+2Q)+cos(4ri«-2a>-2Q)} 

+}|i^-H«V'-|F>^+Tk^'ncos(3ri^-5fi>+2Q)+cos(5n^- 



- (i^ - i^ - -^eV*} {cos {Znt + oi - 4Q) + cos (r^e + Zw - 4Q) 
+ {A/-A/-lft^}cos(2n^ + 2a>~4gS) 

- t^eV* {cos (6n^ - 2a> - 4^) + cos {2nt - 6fi> + 4Q) } ' 

- iH^V {cos (671^ - 4fi> - 2Q) + cos (471^ - 6a> + 2Q)} 
+ T%e>^ {cos (n^ - 7fi> + 6Q) + cos (^nt - fii - 6Q)} 

- xf^Q^ {00s {nt + 5(0 — 6Q) + cos {bnt + (0 — 6Q)} 

+ -fiieV {cos (^nt - 3a> - 4^) + cos {Znt -1a) + 4Q)} 
+ IH^V {cos (7ri^ - 6fi> - 2Q) + cos {but ~ 7fi> + 2Q)} 
+ T^/ {cos (47ie + 2ai - 6Q) + cos (2ne + 4ai - 6gJ)} 



(92) 
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- {2e» -^^ + ^er*-m^'-i^- W«V'} sin 3(n« - a.) 
+ {f«^ - He* + -A^*} sin 4 (n« - a») + f||€« sin 6 (n< - a») 

- {IF - t'W^' + H«V} sin 5 (n< - 0.) - ilf^^ sin 7 (n« - «) 



+ { ley - |eV^ - {4|eV + ^/} {sin 2 (a» - gj) - sin 2 (n< - g3)} 
+ {^fV - ^ey* - ll«y } {sin (3ne + o) - 4g3) - sin (n< + 3«> - 4gJ)} 

- { A«y* - A^y* - WY] {sin (n« - 6a> + 4Q) + sin (5n« - a> - 4^)} 
+ {¥«V - H«V* - M«V + T^A^} sin (3n« - «» -2g^) 

+ {Vey - HeV - MI«V + i^ey*} sin (n« - 3«» + 2Q) 

- {lieV - IW - W«V + ^"j^} sin (5n< - 3«, - 2^) 

- {H«V - M«V - WeV + ^^y*} sin (3n< - 6«> + 2g3) 
+ {|€V - |«V - HieV + ^f} sin (4n< - 2«, - 2S3) 

+ {fcV - |e^' - Mi«V + ^7"} sin (2n< - 4« + 2Q) 
+-i^ey {sin (2n« - 6a) + 4^) + sin {6nt - 2<o - 4Q)} 
+IH«V {sin (4»« - 6a» + 2Q) + sin (6n« - 4tt» - 2Q)} 
+ A«V* {sin 4 (<o - a) - Bin 4 («« - Q) } 

— yj^y* {sin 6 («> — ^) — sin 6 («< — Q) } 

+ ^f {sin (27i« + 4«> - 6g3) - sin (int + 2o> - 6g3)} 

- ^ei* {sin (7to« -Zw- 4^) + sin (3n« - 7<u + 4^) } 

- Mley {sin (5n< - 7o) + 2g^) + sin (7n« - 6a» - 2Q) } 
+ xl^ff^ {sin (6n< + w — 6S3) — sin (jit + 5tt» — 6S3) } 

- -jf^ {sin (7n« - <a - 6^2) + sin (n« - 7a> + 6gJ) } 



. (93) 
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- {6e» + ^yy<^ + ff^y} COS (n< - a») 

- {8e* - fjef + fey } cos 2 (n< - «») 
4- {9e» +f|e^ + ffeV} cos 3 (n« - a») 

+ {2e* — ^ + |cV} cos 4 (n« - o)) + -1^ cos 6 (n« - <m) 
_ {3e» - .^y:*' + l^eV} cos 5 (»i« - <o) - J^'^ cos 7 (n« - a») 

- {6eV - 3^' + fJeV + ^fV} cos (n< + a) - 2Q) 

-i- {4eV - 2ffy - J^eV + ^^} cos (3n< - fl» - 2Q) 
+ {4eV - 2^* - J^i^ + ^^} cos (n« - 3a» + 2a) 

- {«V - ieV - W«y + A«/} cos (3ni - 6«> + 2Q) 

- {«V - i«y - W«V + A^y*} cos (6»« - 3<w - 2g5) 
+ fe*)'^ cos {2ni + 2tt> -4^) - ^^ cos (3n< + 3tt> -6Q) 
— ^e/ {cos (7n« — a> — 6ft) + cos (»i< — 7a> + 6ft)} 

— J^Ley {cos (5n< - 7tt» + 2ft) + cos {Irtt -5a>- 2ft) } 
- 1^* {cos (3nt + a» - 4ft) + cos (nt + So)- 4ft)} 
— 1^* {cos (Znt - 7a» + 4ft) + cos (7ni - 3a> - 4ft)} 
+ AA* {cos (2n« - 6(u + 4ft) + cos (6n« - 2o> — 4ft)} 
— ^V {cos (2n« - 4o) + 2ft) + cos (4»« - 2a» - 2ft)} 
+ J^V {cos (4n< - 6a> + 2ft) + cos {6nt - 4o> - 2ft) } 
+ W«V {cos {nt - 5<w + 4ft) + cos {5nt- a> - 4ft)} 
+ -^ge}* {cos (ni + 5«> — 6ft) + cos {5nt + a> — 6ft) } 

- }eV* {cos 4 («> - ft) + cos 4 (n« - ft) } 
+¥«V {cos 2 (tt» - ft) + cos 2{nt- ft)} 



(94) 



DEVELOPMENT OF THE ELLIPTICAL MOTION. 
f{ntf = {20e* + ^^^ + ^^^] sin {nt - oi) - J^ sin 2 (nt - w) 

-{10^-^e^ + ^^'} Qm3{nt-a)) + 15^cos4:{nt-a)) 

+ ^e^ sin 7 (nt -(o)+ ^eV {sin 2 (a» - ^) - sin 2{nt-Q)} 
+ ^eV {sin {2nt - 4fi» + 2Q) + sin {Ant -2(0- 2Q)} 

I 

~ JeV {sin (4n^ - 6fi> + 2Q) + sin {6nt - 4fi> - 2Q) } 
+ ^eV {sin (5n^ - 7ft> + 2g3) + sin (7n^ - 5fi> - 2Q)} 
+ i^eV^ {sin {nt - 3fi> + 2^) + sin {3nt - a> - 2Q) } 
-■ Wy^ {sin (3n« - 5fi> + 2Q) + sin (5n^--3a>-2Q)} 
+ ^^* {sin (3n^ + fii - 4Q) - sin {nt + 3fi> - 4Q) } 
— ff^V* {sin (n^ - 5a> + 4^3) + sin {5nt - o> -4Q)} 
+ ^eY {sin (3n^ - 7fi> + 4^) + sin {7nt ~ 3fi> - 4Q)} 
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(96) 



/(n^)« = 20e« - ^e^ cos (n^ - ft>) - 30^ cos 2 (n^ - w) 

+ ^e! cos 3 {nt — o>) + 12e* cos 4 (?2^ — co) — 4^e^ cos S (ti^ — a>) 

— 2e^ cos 6 (n« - tt>) + f e^ cos 7(n^ - ft)) - 30eV cos(n^ + tti - 2g3) 
+ ^eV {cos (3n^ - fi> -- 2SJ) + cos (rz^ - 3a> + 2Q)} 

— 9eV^ {cos (5ne - 3tti - 2Q) + cos {Znt-bw + 2Q)} 
+ |eV^ {cos (77i« - 5ai - 2Q) + cos (Sne - 7fi> + 2Q) } 



>■ • 



(96) 



/(w^)^ = 70e^ sin {nt - a>) - 426^ sin 3 {nt - ft)) 
+ 14e* sin 6 (rii ~ a>) — 2e^ sin 7 (n^ — fl> ) 



(97) 
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Ttdt 

{|e» + ie» + ^lyey* + T^^ - T^* - T%eV*} sin (n« - o.) 

- {!«* - A** + tW - 1!^ - A^y* - M«V1 sin 3 (7i< - 0.) 
+ {l^e* + 1^ + ^^«} sin 4 (n< - a») + fiff e« sin 6 (to< - fl*) 

- {Me" + !«' + HI«V'} sin 5 («< - «») - Hf e^ sin 7 (ni - a.) 

- {H«y - 4ieV* - MeV + T^y'} sin 2 (n< - Q) 

+ {i^y' - i^y* - T^«^'} sin 4 (n« - ft) + ^/ sin 6 (ne - ft) 

+ {^- i€y* + 3^;iV«/ - tW- A«V + 1^*} sin(n« + «» - 2ft) 

- {i«y* - i«y* - A«y + A«V* - Trfir«y + lire}'*} sin (n« - 3tt» + 2ft) 

- {|ey« - fey* - |eV + |€»y* - *eV + 1¥5«)^} sin (StK - fi) -2ft) 

+ {Wf- ii«V - T^eV + ihf} sin (2n« - 4«, + 2ft) 
+ { V eV - H«V - |f«V + K^/} sin (4n« - 2o> - 2ft) 

- {M«V - i^«V* - m<^ + T5Vir«/} sin (3n< - 5«> + 2ft) 

- { W «y - iM«V - M«y + ift^«y*} sin (5n< - 3^ - 2ft) 

- {^/-iVy*-A«V}8in(2»«+2a>-4ft) +^jf^«sin(3n< + 3tt»-6ft) 

- {TVey* - tW - Tft^'} sin (n« - 5tt» + 4ft) 
~ { A«y* ~ A«/ ~ IH«y } sin (5n< - <w — 4ft) 
+ { Aey* - tW - l^«V} sin (3n« + a. - 4ft) 
+ {t^* - iW - T^Cy'} sin (w« + 3«> - 4ft) 

+ Hj^* sin (6n« - 2o» - 4ft) + t^^eV sin {2nt - 6«> + 4ft) 
-iH f«V' sin (7n« - 3a. - 4ft) - H^eV* sin (3n< - ?«> + 4ft) 

- tHW sin (n« - 7ft» + 6ft) - ^^e/ sin (7n« - a* - 6ft) 
+ Th^y* sin {Tvt + 6a> — 6ft) + ^^* sin (5n< -Ho* - 6ft) 

- IMK/ sin (7»« - 6<w - 2ft) - ^f^ sin (5n< - 7tt> + 2ft) 

- -gij/ sin (4»« + 2tt> - 6ft) - y^/ sin (2n« + 4o> - 6ft) 

+ fli«V8in(6n<-4a>-2ft) + f|Jey8in(4n<-6«> + 2ft) 



>■> 



(98) 
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^d^finty _ 



v?de 

- {e» + H€» + ^/IjV + 1^«' -■ ifW - MeVI sin (n< - «») 
+ {3e^ + 1^ + J^*} sin 2 {nt - to) 

+ {3e» - We* + A«y* - H-U e' - A«/ - !H«V*} sin 3 (n< - w) 

- {6e* - ie« + HeV} sin 4 (ni - fl*) - Vy^ e» sin 6 (n« - tt») 

+ {W^e*-MM«' + V)ft?«^*}sin5(n«-a,) +ffjf^8in7 («<-«») 

+ {|ey - ieV - iif«V + tI?/} sin 2 (n« - g^) 

+fje^*sin 4 (n« - ft) - T^y* sin 6(n« - Q) - ^^^^sin (2n< -6o> + 4Q) 

- {!!«)'* - ll«y* - T^^'} sin(3n« + «» - 4ft) - fff e*)'' sin(6n< - 2a»-4ft) 

+ A^y* - A^y* - ift/^*} sin (^^ + 3<o - 4ft) 

+ TSW - TJ^ - Tfife«y } sin K - 6a» + 4ft) 
+ {M-T* - M^y* - '^eV*} sin {5nt - a» - 4ft) 

- { AeV - T^«V* - tWW + Tfir^y*} sin {nt - 3^* + 2ft) |- ' ^^^^ 

- {ffeV - ||«y - fHeV + ^ey«} sin (3n< - a» - 2ft) 
+ { W«V - W«y - 'VWeV + W^y*} sin (5n« - 3«> - 2ft) 
+ ||«V - IFy* - 4M«V + TsSW} sin (3ft« - 6«» + 2ft) 

- {2eV - eV* - W«V + -5^/} sin (4ni - 2a> - 2ft) 

- {ieV - ieV - HieV + Tk/} sin (2n< - 4a» + 2ft) 

- ^eV sin (4n< - 6«> + 2ft) - SW«V sin (6n< - 4o» - 2ft) 

- T^/ sin {2nt + 4.(0- 6ft) + ^f sin (4to« + 2a» - 6ft) 

+ fll«*y^ sin {7nt - 3tt> - 4ft) + l^eV sin (3n« - 7<» + 4ft) 
+ fl^feV sin {5nt -7<o + 2ft) + ^^^^ eV sin (7n< - 6a> - 2ft) 
— ^y^ey*sin(5n< + a» — 6 ft) + ^^^* sin (n< + 5a> — 6ft) 
+ MW sin (7»« - «) - 6ft) + yf^ey* sin (n< - 7<u + 6ft) 



24 



THEORY OF THE MOON'S MOTION. 



Th Cvv 

+ { W + IH«^ + Mi«V'} sin 3 (n« - a>) 
+ {^« - J^e« + e^*} sin 4 (n< - o») 

_ {l^V - \%%\^<? + ^W^*} sin 5 i^ - "») 
+ 1|^ Bin 6 (n< - w) - -i^^^ sin 7 (n« - w) 

- UeV - -JeV + ^jVjeV + tI^*} sin (n« + 0, - 2a) 
+ {f <y - ieV - -Jj^iyL^ + ^yV«/} sin (3n« - w - 2^) 
+ {i«V- A«y - im«V + xk^y*} sin (n« - 3a, + 2Q) 

- {|eV - A«V - ^i^VV^ + T^ey*} sin (3n< - 5«; + 2a) 

_ { Jj^^^ _ J^^* _ i^||ieV + Hfey*} sin (5n< - 3tt, - 2Q]! 
+ feV* sin (2n« + 2a, - 4Q) - ^^V«y* sin (3n< + 3a,- 6^2) 

- ^<^ey» sin (7ni - a, - 6Q) — i-sVbV sin (7i< - 7a, + 6Q) 

_ JjiyjjjJ^eV sin (5n« - 7a, + 2a) - ^^^^yi^^ sin (7n< - 5a» - 2Q) 

- fffeV sin (3n« + a, - 4a) - tVW^ sin (n< + 3a, - 4a) 

- fff«V sin (3n< - 7a, + 4^) - ^VW«V sin (7«< - 3a, - 4^) 
+ T^yeVsin {2nt - 6a, + 4^) + f|«V sin {6nt - 2a, - 4^) 

— "j^ey sin (2n« - 4a, + 2^) - 17eV sin (4n« - 2a, - 2^) 
+ V«V sin (4n« - 6a, + 2^) + ^^y sin (6n« -4a, - 2^) 
+ tA^I^V sin (n< - 5a, + 4a) + ^^y^gy sin (6n« - a, - 4Q) 

- 2^* sin 4 (»< - a) + H«V sin 2 (n< - a) 

+ sh^ sin (nt + 5a, - 6a) + HIV sin {5nt+ a, - 6a) 



(100) 
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tk^^S- = {i^ + t!I7«' + i/V«V1 sin (n« - «.) - fe* sin 2 (n< - «,) 
nav 

- {W - ¥i^' + W«V*} sin 3 {nt - a») + 32e» sin 4 (n« - a») 

+ { W^* - ^*W¥^^ + IM«V'} sin 5 (n« - a») - V e* sin 6 (n< - o.) 

+ HfM^^ sin 7 (tK - a>) - |ey sin 2 (n< - g3) 

+ feVsin (2ni - 4a> + 2Q) + V^V sin (4n< - 2tt» - 2gJ) 

— I^V sin (4n« - 6a» + 2Q) - A^eV sin (6n< - 4tt» - 2Q) 

+ W^«V sin (6n< - 7«» + 2gJ) + J-ff^tgy sin (7n< - 5a> - 2^) 

+ t^«V sin (n« - 3tt» + 2a) + ^VW^* sin (3n< - « - 2Q) 

- W«y sin (3ne - 5o. + 2^) - ^H^eV sin (5n< -Zat- 2£3) 
+ W«V* sin {Znt + o» -4^) - xfr^V* sin (n« + 3tt» - 4g2) 

— fffjeV^ sin (n< - 6m + 4g3) - ^^e^* sin (6n« - w - 4gJ) 
+ ^^^ey sin (3n< - 7«» + 4S^) + ^^H^ sin (7n< - 3a> - 4Q) 



^; 



(101) 



7 ^0 \j1 ~ "^ "A^^ ®^^ (^^ — co) + fe* sin 2 (n^ — w) 

- ^W«^ sin 3 (n<- a») - ^e» sin 4 {rd - a») 
+ 4|ic^ sin 6 (n< - w) + ifie* sin 6 (n< - w) 

+ ,ifj«V sin (ni + a> - 2Q) - -^eV sin (3n< - fl» - 2g3) 

- .^^ sin (7i« - 3tt» + 2Q) + ^eV sin (5n« - 3«> - 2^) 

+ W«V sin (3n< - 5a> + 2Q) - J^lfVy sin {Int -5a>- 2Si) 

- W«V sin {5rd -7to + 2Q) - J-ff^e'^sin 7 (n« - a») 



>> 



(102) 



- A|^^ sin 5 (n< - a») + ifH^^ sin 7 (n« - o») 



/• 



(103) 



If we now substitute equation (91) and equations (98-103) in equation (90), 
we shall obtain the following value of the true longitude v of the moon in terms 
of the mean longitude nt. 

6 
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v=nt 

+ {2e-y' + ^-iey*+\ef + i^' + :^*} sin {nt-a>) 
+ {f «* - tt«* + i^y* + 1^ - i«y -1^/} sin 2 (n< - «,) 

+ { W«' - ififi* + H«V*} sin 4 (ni - a») + VWe* sin 6 (rK - o») 
+ { VWe» - fMF + itie^*} sin 5 (««-«,) + efH«' sin 7 (n« - 0,) 

- {if-\y*-<^ + yy +neY + ^/} sin 2{nt-Q,) 

+ {A/ - isf -Wr*} sin 4 (nt -Q)- ^/ sin 6{nt-Q) 

■ + (ix* - iy* - i^V^V - li^V + A/ + T^«V} sin 2 («, - a) 

+ (ifV/ - A/ - AeV*} sin 4 («» - a + y^/ sin 6 (a; - gj) 

+ {iey» - iey* - T^ijeV + ^^« + iey* - T^l^^} sin (n« + a» - 2a) 

- {iey» - i€y* - H«V + HeV* + tW - Hi«V} sin (3n< - fi, - 2^) 

- {H«y -il«V - W«V + ^f} sin (4n« - 2fl* - 2£3) 

- { Ary* - T^/ - H«V} sin (2n« + 2io- AQ) 

+ { A«V -• 5')r«V + T%«y } sin (n< - 3a» + 2gJ) 
-{Jey*-iV— |«V} sin (3n« + a» -4^3) 
+ {-i^y* -W~ U^*} sin (»i« + 3«» - 4a) 
+ {ie/ - i«/ - W«V} sin (5n< - a» - 4^) 

- {||«y - M«V + M* + Wi?eV} sin (5n« - 3a» - 2Q) 
+ (k^V sin (2n< - 4a> + 2^) + l^eV sin (6n« - 2«» - 4Q) 

- W «V sin (6w< - 4a> - 2a) + ,^/ sin {4nt + 2«> - 6^) 

eV sin (3n< — 5a> + 2a) 



- -^f sin (2n< + 4a> - 6a) - 

— i^ sin (6n< + 0* - 6a) + H«y* sin (7i< + 6tt» - 6a) 
— ■^' sin {7nt — « - 6a) - -^e/ sin {Bnt + 3«> - 6a) 
+ {Uei* sin {7nt -3a>- 4a) + ^Wi^^ sin (7n< -5«> - 2a) 
+ W«V sin {5nt - 7«» + 2a) + tI^* sin {nt-5(o+ 4a> 



(104) 
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6. Having found the value of v in terms of nt, if we now substitute it in equa- 
tions (73) and (39), we shall also obtain the values of r and in terms of nt For 
this purpose we shall first develop equation (73) in a series, carrying the approxi- 
mation to terms of the fifth order, and we shall find 



T 

+ ri^* cos (v — 6ft> + 4S3) + yfr^^ cos (pv—o) — 4S3) 

— yf^^cos(3i; + w — 4Q) — yf^ey*cos(v + 3a> — 4SJ)+ ^ey*cos 3 (v— o)) 



\ . (105) 



In order to substitute the value of v in this equation, we shall observe that if 
we put all the terms of the second member of equation (104), except the first, 
equal to j9, we shall have 

v = n^ + j8. (106) 



This gives 



mv + a=^ (mnt + a) + mj8. 



(107) 



m being any positive whole number, and a any angle whatever. 
Equation (107) gives 

cos {mv + a) = cos (mnt + a) cos m^ — sin (mnt + a) sin mj8. 



(108) 



Since ^ is a small quantity, we may develop sin mj8 and cos m^ in a series, 
and we shall have, with sufficient accuracy, 



sin m^ = m^ — -Jm^j^, 



cosmj8 = l-^m"^ + i^*^. 



(109) 



We shall then find, as far as terms of the fourth order, 

^=2^ + '^^ + ^f + ^^coa{nt-(o) 

— { 2^2 - |e* + ^/} cos 2(nt-(o)—^ cos 3 (nt — oi) 

-.^e* cos 4 (n^ - o)) - fJeV cos 2 (n^ - Q) + fJ-eV cos 2 (fti - Q) 

— -^y^coa^nt - Q) - ^y^cos4(fti - Q) +3^y*cos (2n^ + 2fti -4g3) 

— ef co8(nt + a)-2Q>)+^ coa{3nt -w-2Q,) + ^ cos {nt-3a)-^2Q) 
+ 1^ cos (471^ - 2fti - 2Q) + 3^6^ cos (2n^ - 4a; + 2gJ) 



(110) 
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^ = Ge* sin (nt — tt») +^^Bin2(nt — to) — 2^ Bin S{nt — at) 
- ^ sin 4 (n< - to) - feV sin 2 (n« - gj) + |eV sin 2 (fi» - gj) 
+ |ey sin (4n< - 2«» - 2£3) + |«y sin (2n< - 4<u + 2a) 



ail) 



^ = 6e* — 8e* cos 2 {nt — a>) + 2e* cos 4 {nt — a>) . 



(112) 



Therefore equation (108) will become 

cos {mv + a) = cos {mnt + a) X 

{l-m"(e" + ^ + ^y*)+imV-|mVcos(7i«-a>) 

+ JmV cos 3 (nt — 0)) + (ff^V +t*5^V) cos 4 (n^ — w) 
+ l^m'eV cos 2 (n< - gj) - iim^eV cos 2 (a> - gj) 
+ ^mV cos 4 (n^ — Q) + -gfym V cos 4 (ft> — ^) 

— -gV^V cos (2n^ + 2a> — 4Q) + ^efm? cos (n< + co — 2S2) 
— Jm^ey^ cos (3n^ — a> — • 2Q) — ^^ef cos (n< — 3ft> + 2Q) 

— fJm^eV cos (47i« - 2ft> - 2£^) - -^^e^ cos (2n^ - 4ft> + 2gJ) } 

— sin (mn< + a) X 

{ (2me- Jme^-77iV)sin(n^— o)) +(}m6^---J^me*+3^my*--|mV)sin 2(7i^— co) 
+ (Ifme^ + "ImV) sin 3(n^ — «>) + {^^^me^ + -frnV) sin 4 (n< — a>) 

— (Jmy^ ~ -^mj^—meV ~|m^6V) sin 2(n^— Q)— J^ey^sin(3?i^— fi> --2g^) 
+ (imy2-|m/~3^ey-|m3eV^sin2(ft>-^)+imej^sin(n^ 

+ ^my^sin 4(7i< — ^) + ■^my*sin4(a> — ^) — ^^my* sin (2n^ + 2fti — 4Q) 
~ (If meV + -Im^eV^ sin (4n< - 2ai - 2a) - •tm'eV sin (2n^ - 4ai + 2gJ) } 



. (113) 



If in this equation we put in succession m = 1, 
a = ft>—2gj; m = l, a = — 3a) + 2£J; and m = 3, a = 
get the following equations : 



a = — a); m=l, 
tt> — 2SS; we shall 
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+ {e- fe*} cos 2 {nt - co) + {fe^ - ff|e^ + fy^} cos 3{nt - w) 
+ fences 4 (?^^— o)) + f|fe* cos 6 (n^ — oi) 
~{iy'-i^/-M^}cos(3n^-a>-2Q) 
+ {iy" - i/ - A«^'} cos (n^ + «> - 2gJ) 

+ xf^y^cos {5nt — w — 4SJ) — y^y* cos (Sn^ + fi> — 4SJ) 
+ xfry* cos (n^ + 3fti — 4Q) — yj^y* cos (n^ — 5(0 + 4S3) 
+ -fey^ cos 2 (n^ — Q) — fey^ cos {4:nt — 2a> — 2Q) 
- ley^cos {2nt - 4a> + 2Q) - ^ey^ cos 2{a)-Q) 
-fi^Vcos (5nt - 3a> - 2gJ) - -^eVcos (3n^ - 5a) + 2^) 



; (114) 



C0B(i; + fii - 2gJ) = (1 -6^ cos(n« + o) -2gS) + fe^ cos(37i^ - iw - 2g^) 

- ^ cos ('/i^ - 3fti + 2£^) - ly^cos (3/1^ + a> - 4g3) 
+ ^y^ cos (n^ + 3^ — 4g3) + e cos 2 (n^ — Q) — e cos 2(ft>— Q) 



; (115) 



cos(v -3ai + 2^) = (1 -e^ cos {nt-3(o + 2S3) -\e^ cos (n^ + ft>'-2gj) 

+ |e* cos (3n^ - 5a; + 2Q) + ^f cos (n^ + 3ft> - 4gJ) 
■~iy* COS (n^ — 5fi> + 4S2) + e cos (2n^ ~ 4a> + 2Q) 



I ; (116) 



— ecos 2(a> — Q) + •Jy'cos (nt — o)) — -^y^cos 3 {nt — ft>) 



cos(3v-fi>-2Q) = (l~9e^cos(3n^-a;-2Q) +-^e2^og(5^^_3^_2ga) ^ 

+ ^e" cos (?i« + fli - 2S3) + f / cos (3n^ + ft> - 4g3) 

— f y* cos (5n^ — oi — 4gJ) + 3e cos (4n< — 2ft> — 2Q) 

— 3e cos 2 (nt — ^) + f y^ cos (n^ — a>) — |y* cos 3 {nt — (o) ^ 



. (117) 



If we now substitute these values in equation (105), and at the same time 
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change v into nt in the terms of the fifth order in that equation, it will 
become 



^ — ^ = l-e2 + e{l-|«" + tVV«*}cos(n«-a)) 
r 

+ e" {l-^e"} coQ2{nt- w) + {l^ -mef"} coBZ{nt-a)) 

+ f e* cos 4 (nt — a))+ f|^e* cos 6 {nt — (o) 



(118) 



In the substitution in equati6n (105), the coefficients of all the terms contain- 
ing Q, become identically equal to nothing, and the quantity y disappears from 
the coefficients of all the remaining terms; whence it follows that the radius 
vector of the orbit is entirely independent of the position of the orbit. We shall 
therefore have 

^ = l + e{l-ie' + j^e^}cos(nt-a})+e'{l-'ie'}coQ2{nt-a}) 1 

+ {|e^--^e'} cos 3 (nt-o)) +|€*co8 4 (n«- a>) + fUe'cos5(ne- co) j 

In order to find the tangent of the latitude it is only necessary to find 
sin {v — Q), and multiply it by y. To find sin (v — Q) we shall observe that if we 
change cos {mnt + a) into sin (mnt + a), and — 8in(mn^ + a) into + cos {mnt + a), 
cos {mv + a) will change sin {mv + a). If we then make these changes in equa^ 
tion (113) and put m = l, a = — ^, we shall get, after multiplying by y, the 
following value of y sin {v — Q,), or tan d, 

-{if''Thf-U^}^^^^{nt-Q)+j^y'sm5(nt-'Q>) 

+ 6y { 1 — fe^ — If} sin {2nt — a} — Q)—ey sin (a) — Q) 

+ {|^ -■ ^V + Ay* - A^} sin (3n^ - 2ai - gj) 

+ {i^-if-ii^'Y + 'iif + -heY}sm(nt-2o, + Q) 

+ {iy'-A/-H^}sin(n^ + 2a>~3g3)-i6r^sin3(ai-gS) [.(120) 

+ f 6^ sin {Ant — 3a> — ^) + f e^ sin*(4n^ —Scj — Q) 

+ {-jV^ -W} sin {^nt -3(0 + 0;)+ fff eV sin (5n^-4a; - Q) 

+ {yf^V - A^} sin (3n« - 4fi) + Q) + yf^y* sin (n^ + 4a> -5Q) 

+ { A^ - ihf] sin {nt - 4ft> + 3S3) - ^ef sin {Ant -a)- 30) 

-^y«sin(3n^+2ft>-5gJ)+^sin(2n^+co--3gJ)"|Jeysin(6n<--2a>-3a) J 
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This equation will give 

+ -^f Bin 5 (nt- Si) -^ef Bin (a)- Q) + ^ Bin {2nt + a» -Sgi) 

+ Je}^Bin(2n« - o» - Q)- Jey*8in(4TO<- a> - 3Q)- ^y»8in(3n<+2tt»-5g3) 

+ i^f{f + 9^} sin (Snt - 2a» - gS) - ^^eV sin {6nt - 2tt» - 3Q) 

+ -hf{^-f}Binint-2a> + a) + -^f{f-7.^}Bin{nt + 2w-3a) 



^ . (121) 



ita.n'^e=\f Bin {nt- Si)- ^/BinB {nt-Q)+-^f Bin 5 (nt-Q). 



(122) 



Now we have 



^ = tan (? -Jtan'tf + ^ tan**^ ~, etc. 



(123) 



If we now substitute these values of tan d and its powers in equation (123) we 
shall get the following expression for the latitude of the moon, which is correct to 
terms of the fifth order : 



d=r{i-^-if+i^* + -^y*+U^}Bin{ru-gi) 

- {-hf - Thf - Wf} sin Hnt -0.)+ ^f sin 5 (nt - Si) 
+ ey{l-ie'-ff}Bin{2nt-(o-Q)-ey{l-if}Bin{a>-a) 

+ {|«V - H«V + At* - H«y } sin (3n« - 2«» - g3) - iey» sin 3 («> - gj) 

+ ii^ - if - is^y + ii/ + -iieY} Bin {nt- 2a, + Q) 

+ {if -iff- if «V} sin (nt + 2a,- SSi) + |cV sin (Ant - 3«» - Q) 

- ief sin (Ant - o» - 3^) + ^ sin (2nt + w — 3Si) 

+ {^JjeV - ief} sin (2nt - 3a> + Q) + Hf eV sin (5nt -Aw -Si) 
+ iih^ - A«V} sin (3n« -Aa, + Si)+ jhsf sin (nt + Aa}-5Si) 
+ { A«V - ihf} sin (nt -Aa, + ZSi)-itf sin (3n< + 2<o -5Si) 

- HeVsin (5nt -2a>- SSi) 



(124) 



7. We have thus found the values of the three corordinates r, v and d in terms 
of the time. We may however find d dir^ectly from the differential equation (28), 
which will serve as a verification of all the developments which have thus far 
been made in the determination of t, v, r and d. For this purpose we shall sub- 
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stitute the values of c, c' and c" in equation (28), by which means it will 
become 






(125) 



Now equation (119) will give 



a 



^ = l+-^ + |e*+{2e + |e» + |4e*}cos(n«-ft>) 
+ {fe^ + |e^} cos 2 (rwt - ft)) + J^* coa 4 (n^ - fi>) 
+ { W - M^} cos 3 (n^ - fl>) + J^^fe* cos5 (n^ - fl>) J 



>i 



(126) 



and equation (113) gives, by putting m = 1, and a = — S3, 



C08(»- a) = {l-^ + ^f ■^■^-■hf-'h^} cos(n«- G) 

— ecoa{(n — Q) + {e — ^ + ^«y*} cos {2nt — m — Si) 

+ {l«* - tt«* + A/ + A«V} cos (3n< - 2a> - d) 

- (i)^ - rky* - li«V} COS 3 (n< - gj) 

+ {if—^y*-U^}coa{nt + 2a>-SSi) 

+ f e» cos (4n« - 3a> - g3) - {^ + -J^} cos (2n< - 3<u + a) 

— -J^y* cos 3 (tt> — Q) + ^ey* cos (2rrf + a» — 3Q) 

— ley* cos (4n< - a> - 3Q) + ff|«* cos (5«< — 4oi — gj) 

— {-A^* + A«V} cos (3n«— 4a> + Q) + xl7y*cos 6 (rrf— a) 

+ ihy* cos (w« + 4«> - 5Q) - {j^y* + ^^jeV } cos (nt - 4tt» + 3g3) 

- -^f cos (3n« + 2a> - 5Q) - 1^ cos {5nt -2<o- ZQ) 



►. (127) 



If we now multiply equations (126) and (127) together, and the product by y, 
we shall obtain 
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- (i/ - T^y* - M«V } cos 3 (n< - S3) + T^Tr/ cos 5 (n< - a) 

+ {2ej' - |eV + i*)^} cos (2n« - tt» - Q) + J^ cos (4n< - 3a» - a) 

+ {¥«V -■ ¥«V + H«V + Ax*} cos (3»< - 2a> - Q) 

+ (i/ - Ax* - ii«V} cos (»< + 2fli - 3£J) 

+ {i«V - i)^ + A«V + Ay* + A«V} cos (Ji< ^ 2a> + Q) 

- Uey* - i«V} cos (2n« - 3a> + Q) + ^ cos (2ni + a» - 3gJ) 

+ {^!^V - H«V} cos (3nt - 4oi + gj) — J^cos (4n< - a> - 3g3) 
+ {Vt«V -^y*} cos(n< -4«» + 3Q) + ^Vy^Vcos (5?i< -4fi> - a) 
+ Tfyy* cos (nt + 4a> - 5a) - -^Z cos (3n< + 2a» - 5gJ) 
— IJcV cos (5n< - 2a» - 3Q) 
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. (128) 



If we multiply this equation by 

v/l + y* 



and put ^-^ = n, we shall obtain 
a* 






= y { 1 - e* - fy* + 1^/ + t/V«* + H«V} cos (n< - g3) 

- iy* {1 - Hy* - ¥«*} cos 3 (n< - ft) + tI^T* cos 5 (»< - SS) 

+ 2c)' { 1 - |e* - f y*} cos (2n< - a> - gj) + J^cV cos (4n< - 3tt» - a) 

+ {¥^ - Wy - ■|i«V + A/} cos (3»i< - 2© - g^) 

+ (iy* - Ay* - if«y } cos(n< + 2a> -3Q) - ||ey cos(5n« -2a> -3g3) 

+ (Fy - i^* - A«V + liVy* + 7^V«V} cos (n< - 2o* + g2) 

— {i«y* - ie^} cos (2n< - 3a» + a) + ^ cos (2n< + tt» - 3a) 

+ {t^V - H^V } cos (3n« - 4a> + a) — ley* cos (4%< - a> - 3Q) 

+ {^cV -yky*} CCS (w< -4«» + 3a) + ^^Vcos (5n<- 4a> - £1) 

+ ify/ cos {nt + 4tt» - 6a) - -A/ cos {Znt + 2a> - Sa) 
7 



(129) 
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Equation (129) gives, by integration, 

(? = ((?) +y {1 - 6« - |y« + T^« + T^/ + Hey} sin (n< - g5) 

- Hif-Thf - iW} sin 3 (nt -0,)+ jihif sin 5 (nt - Q) 
+ ey {1 - Je" - 1/} sin (2n< - «> - Q) - ^|ey sin (5n< - 2tt» - 3^) 
+ {i^ - H^ + ^/-UeY} sin (Znt- 2a, -O) 
+ {\e^-\f-ise*Y + iif + -h^}sin{nt-2a, + Q) 

+ {if --hf- Wf) sin (^< + 2a> - 3Q) + f eV sin (4n< - 3a> - gj) 
- -Jcy* sin (4n« - «tf - 3g3) +i«y»sin (27i< + «> - 3Q) 
+ {^«V - i*?^} sin (2n« - 3«* + Q) + flJeV sin (5n< - 4a» - gj) 
+ {xb^V - A«V} sin (3n< - 4a< + gj) + yf^y» sin (n< + 4a> - 5S3) 
+ {-Mf - ikX*} sin {nb - 4a> + 3Q) - -g^/ sin {Znb + 2«» - SQ) 



> • 
> 



(130) 



(^) being the constant quantity to complete the integral. It is evident that the 
two expressions for the value of given in equations (124) and (130) will be iden- 
tical if we make the constant 



{0)==-ey{l-\f}^m(co-Q^)-\efB^nZ{(o-a, 



(131) 



which makes equation (130) satisfy the condition that the latitude of the moon 
shall be equal to the latitude of the perigee of the orbit when nt = Wy or when 
the moon is at the extremities of the transverse axis of the orbit. The perfect 
agreement of these two determinations of the value of proves conclusively that 
all the preceding analytical developments have been correctly made. 

8. In order to show, by a few numerical examples, that the values of v, r and 
Oj given by equations (104), (119) and (124), are correct, and at the same time 
show that the development of these quantities in series to terms of the fifth order 
is sujficient, in the lunar theory, we shall now reduce these equations to numbers 
by using the values of e and y, corresponding to the elements of the moon's orbit 
which were employed by Delaunay in reducing his equations of the lunar theory 
to numbers. These values are, e = 0.05489930, and y = 0.09004560. And 
instead of giving the value of r in terms of the moon's mean distance as the unit, 
we have multiplied equation (119) by the constant term of the moon's parallax, 
supposing it to be equal to 3422."3. 
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We shall therefore find 



v=nt + 22638."97 sin (n« - w) + 777."074 sin 2 {nt - a») 

+ 36."997 8in3{nt-a}) +2."010sin4:{nt-(o) + 0."118ain5{nt-a>) 

- 411."374 sin 2 {nt -Q,) + 415/'469 sin 2 (a> - Q) 

- 45."255 sin {Snt - oi - 2Q) + 45."691 sin {nt + m- 2Q) 

- 4."096 sin {Ant - 2a> - 2Q) + 0/'093 sin (n< + 3a> - 4g5) 
- 0."186 sin {3nt + w- 4Q) + 0."093 sin {5nt - o. - 4Q) 

- 0."340 sin {5nt -Sto- 2Q) - 0/'848 sin (2n« + 2«» - 4Q) 
+ 0."424 sin 4 (n< - Q) + 0."424 sin i{(o-gi) 



>• (132) 



= 18461."23 sin {nt -0,) + 1012."716 sin (2ni - a> - Q) 
- 1017."590 sin («>-«) + 62."519 sin {Znt - 2(o - Q) 
+ 18."584 sin {nt + 2(o- 3Q) - ll."662 sin {nt -2a> + Q) 

- 5."993 sin 3 {nt -0) + 2."067 sin {2nt + w-Sgi) 
- 1."0335 sin (4»< -m-SQ,)- 1."0335 sin 3{a}-Q.) 

+ 4."098 sin {Ant - 3a> - Q) + 0/'2746 sin {5nt - 4tt» - ft) 
+ 0."0286 sin {nt +4a» - 5^) - 0."1206 sin {5nt -2a)- 3Q) 

- 0."0191 sin {3nt + 2(0- 5Si) - 0."0024 sin {nt -Am + SSi) 
+ 0."0057 sin 5 {nt -Si)- 0/7773 sin {2nt -Sa> + Q) 

- 0/'0520 sin {3nt -Ata + Q) 



J 



(133) 



It = 3422. "300 + 187. "811 cos {nt - to) + 10."304 cos 2 {nt -^o)^ 



+ 0."636 sin 3{nt-to)+ O."0414co3 4 {nt - w) 
+0."0028 cm 5 {nt-io) 



> ; 



(134) 



TT denoting the moon's parallax. 
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"We shall now compare these approximate values of v, and it with the values 
derived by an exact calculation from rigorous formulae for the same co-ordinates. 
For this purpose we shall observe that if we denote by u the eccentric anomaly, 
and by nt — o)' the mean anomaly measured on the plane of the orbit, the rela- 
tion between u and tiit — iof will be given by the equation 

92^ — co' = t^ — e sin t6. (135) 

Then we shall have the following equation for the determination of the true 
anomaly v' — «>', and radius vector r. 



tan •J-(i;' — a>') = ^ . tan \ti 

r = a(l — e cos v) 



> • 



(136) 



Equations (135) and (136) are entirely rigorous ; and if, for a given value of 
the mean anomaly nt — o)', we compute the corresponding value of u with accu- 
racy by means of equation (135), we can then compute the corresponding values 
of v' — a>' and r by means of equations (136). Now we also have the value of 
ft>' — S3 by means of the equation 

tan (a>' — Q,) = tan {(o — S3) sec i ; (137) 

then adding the value of v' — co' to o)' — Q, we shall have the true distance, 
v' — Q>j of the moon from the node measured on the orbit. 

Then the spherical triangle vv'Q> (see figure, Art. 10, Introduction) will give 



tan (v — S3) = tan {v' — S3) cos i ' 
and sin d = sin (v' — S3) sin i 



(138) 



The formulae (135-138) will give the rigorous values of r, v and d correspond- 
ing to any value of the mean anomaly nt — w'. 

li we suppose that i = 5° 8' 43."28, and ft)-S3 = 30°, we shall find that 
o)' — Q = 30° 6' 0."99. Then if we designate the values of v and 0, derived from 
equation (138), as "observed values," and the values of the same quantities 
derived from equations (132) and (133) as "calculated values," we shall obtain 
the following table showing the values of the co-ordinates v and for the values 
of the mean anomaly given in the first column of the table. 
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Mean 
Anomaly. 


V — 


OfiSE] 


avED. 

d 


Calculated. 


Calculated.- 


—Observed. 


o 

30 


63° 20' 


58.23 


o / // 

+ 4 36 4.30 


63° 20' 58'06 


o 

+ 4 


36 l30 


-d:i7 


o'.'oo 

• 


40 


74 18 


13.61 


4 57 16.05 


74 18 13.72 


4 


57 16.06 


+ 0.11 


+ 0.01 


60 


95 45 


20.32 


5 7 10.39 


95 45 20.51 


6 


7 10.34 


-fO.19 


-0.05 


80 


116 27 


0.21 


4 36 33.14 


116 27 0.22 


4 


36 33.10 


+ 0.01 


-0.04 


100 


136 19 


36.09 


3 33 29.25 


136 19 35.97 


3 


33 29.18 


-0.12 


-0.07 


130 


164 46 


9.63 


4-1 21 18.38 


164 46 9.80 


+ 1 


21 18.40 


+ 0.17 


+0.02 


160 


192 4 


24.48 


-1 4 44.42 


192 4 24.44 


-1 


4 44.40 


+ 0.04 


+ 0.02 



The columns of residuals show that the development of the longitude and lati- 
tude in series as far as terms of the fifth order can in no case be in error to a 
greater extent than about 0.''2 in longitude, and 0."1 in latitude; and it would 
therefore seem that the development of the lunar theory to that degree of approxi- 
mation would be amply sufficient for all the purposes of astronomy, except for 
those terms of perturbations producing the inequalities of long period. We would 
here also observe that the value of the parallax given by equation (134) can in no 
case differ from the rigorous value of that co-ordinate to a greater extent than 
0."1 ; and a more perfect agreement between approximative and rigorous formulae 
could hardly be expected or desired. 

9. The whole of the mathematical theory of the moon's motion is contained in 
the preceding articles (1-8), when we neglect the consideration of the effects 
of the disturbing forces. The expressions for the three co-ordinates v, and r, 
given in equations (104), (124) and (119), are perfectly homogeneous in their 
development, every circular function or argument of the different terms of the 
equations being measured on the fixed plane of projection instead of the plane of 
the orbit. There are therefore no terms depending on the mean distance of the 
moon from the perigee measured on the plane of the orbit, usually called the 
mean anomalies; but the equivalent terms depending on the difference of the 
mean longitudes of the moon and of the perigee measured on the plane of projec- 
tion is given instead. For this reason it is unnecessary to make any reference to 
the distance of the perigee from the node of the orbit measured on the plane of 
the orbit, but simply to designate the longitude of the perigee by the distance of 
its projection from the origin of longitudes and measured on the fixed plane. And 
therefore the expressions for these three co-ordinates which we have given are 
simpler than the equivalent expressions given by Delathstay, who reckons longi- 
tudes on the fixed plane, and the other arguments of his equations upon the plane 
of the orbit, referring them to the fixed plane in terms of certain angular func- 
tions depending on the distance between the node and perigee, and also to the 
mean anomalies which are measured on the plane of the orbit. Delaunay's 
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theory therefore requires the consideration of two planes of reference, one for the 
longitudes and the other for the anomalies. For this reason his formulae are less con- 
veniently adapted for developing the perturbations than the expressions which we. 
have given. We may observe, however, that Delaunay's expressions for v^and d 
Hjay be made identical with equations (104) and (124), as far as terms of the fourth 
order, by changing the y of his notation into ^^(1 — fy^i and also changing g into 
{co — Q)-\-\'f sin 2 (a> — S^). We may then correctly put his Z = ri^ — o), and then 
the two systems of equations will be made identical. Considering the intricacy 
and importance of the subject, and the great extent to which it is necessary to 
carry the approximations, it is very essential to reduce the fundamental elements 
of the problem to as simple a form as possible; and it is believed that the differ- 
ential equations of the moon's co-ordinates, v, and r, given in § 1, leave little to 
be desired so far as the determination of these quantities is concerned. We shall, 
however, in order to render this chapter more complete as a basis for the further 
development of the lunar theory, now^ive the analytical expressions for the vari- 
ations of the elements y, S3, a, e and w^ depending on the action of the disturbing 
forces. 

• 

10. For this purpose we shall observe that we have already found 

tan Q = ^. (139) 

If we suppose that c' and c" are variable, from the action of the disturbing 
forces we shall get, by differentiating equation (139) and dividing by dtj 

do, _ cos^^ fcZc" c" de 






dt c' \ dt c 



But we have c^'^cy cos Q, and c" = cy sin ^ ; therefore equation (140) will 
become 

dSi _^cosQ, dc^' ___ sin Q dcf ,^ ^^ ^ 

dt cy ' dt cy dt* 

If (/ and e" are variable by reason of the disturbing forces, equations (9) and 
(10) will give 

-—- = —(—- 1 tan ^ sin V — (-— ricosv: (142) 

dt \dv) \dd ^ ^ 



de 
"dt 



" (dE\, r. (dR\ . \ ,-. * 

- = — — tan^cosv — (-— Ismv. (143) 

5 \dv \dd) ^ ^ 
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put 



If we now substitute these values of — and — - in equation (141), and also 

at at 



we shall find 



da 



1/TT7 






dt j/a/jL (1 



— — . J sm iv 



— Q>) cos {v — 



('-«)(f)-=^^(f)}- (1**) 



In the same manner we may compute the diflferential variations of the 
elements y, e and w; but the process is simple, and it is unnecessary to give 
it in detail. We shall therefore insert here the equations to be employed 
in the further development of the lunar theory, the solution of which cannot 
fail to give the moon's co-ordinates with all the precision required by obser- 
vation, if the inequalities of her motion are produced by the forces of gravi- 
tation. 

If we therefore substitute the values of c, c', c", / and /', in equations (19), 
(20) and (21), we shall obtain the following differential expressions of the co-ordi- 
nates V, and r : 



dt v^l+j^r^cos^ffl i/a//(l-e2) 



{ 



dv _ \/afi{l-e) (^_ vTT7 



■sm-' 



(B) 






— —^ \ \ \&nO sin v|— — |+cost>(-— - 1 [ dt 

+ -^—— C \ tan/?cosv|— — 1— sint;(-— I \dt 
r* J \ \dv/ \dej) J 



(C) 
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dr 
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dt ^afi{l-^j 



uecoadtcoa06in(v — a>)\l ^ ^ CI — \dty 

\ Vaaii-^J\dv) J 



r VT+7 



+ i 



Vafi (1 - e*) 



-cos 9 cos 



vr|r»{c£»8i 



sin/? sinv — dv cos0 cos 



''(f) 



-I- {2rdv cos /? sin V — 



cosv 
cos^ 



*'(f) 



+ {2rcW cos ^ sin v + dr sin /? sin v} 




_ i/TT? 



l/a/£(l-0 



1^=008/? sin V r J r'fcW sin/? cos v + c?i;co8^sinv}( — ) 
-6^ J 1 ' Adr/ 



+ {2rdlv cos ^ cos v + dr] 

cos^ 

+ {2^d0 costf cosv + dr sin ^ cos v} 








c^t 



; (D) 



—I- = ^ — J sin (v — bd) cos (v — ft) 

dt i/a//(l-6^1 



/c?i?V_ 8in(v 



- ^) /<i. 



d^ 



(E) 



<^y_ i/TT7 



<i< l/a;t(l-e*) 



^{,cc^,»-a)(f)-c»(.-a)(f)|, 



(F) 



/le cos ^0-^;- =r^\ sin ^sin(i; — <o)— — cos d cos (v — (o)— [ ( — — | 
'^ die 1 ^ ^di ^ dt)\dr) 

, fcos(v — ft>) dr o /J • / jiv} ldE\ 

+ \ — ^^ '- — 2r cos ff sm (V — ft>)— - H^r-| 

1 co8» dt ^ ^dt)\dv) 

—i sin^ — + 2rco8& — y sinfv — ft>)(— — | 
1. d^ do ^ \ddj 



(G) 
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fi—=7^\ {sin^ocos^— cos^o8in^co8(v— tt>)} — --cos^ocos^8in(t;— tt>) — [ ( — ) 
at ( at dt) \dr J 



^sin {v — ft>)— + 2r {cos do cos cos(v — a)) + sin Oq sin 0} — I ( — | 

cos^ dt dt J \dv J 



f dv 

{sin Po cos p — cos ^0 sin tf cos (» — «>)} — 

I dt 

— 2r {cos&ocos^cos(v--ft>) 4- sin^osin^ — 

dt 



dR 



dd 



.(H) 



And lastly, equation (18) will give 

We have already integrated equations (B), (C) and (D), when the function JR 
is equal to nothing ; and if we substitute these approximate integrals in equations 
(E-I), we shall obtain the approximate values of the variations of the elements 
of the orbit. "We must then substitute the corrected elements a, e, y, w and Q, 
in equations (B), (0) and (D), by which means we shall obtain corrected values of 
V, d and r, which will contain terms depending on the first power of the disturbing 
force ; then with the corrected values of the co-ordinates v, and r, together with 
the new elements of the orbit, we must repeat the solution of equations (E-I), by 
which means we shall obtain corrected values of the variations of the elements, 
with which we can repeat the solution of the equations which determine the 
co-ordinates v, and r. By continuing the solutions in this manner we may 
obtain all the terms of sensible magnitude arising from the disturbing fojrce, 

11. In bringing to a close this first chapter of a lunar theory, it will be inter- 
esting to compare the results of our analysis' with the conclusions arrived at by 
other mathematicians, in so far as the different methods of investigation will 
permit. For this purpose we shall observe that both La Place and Plana in 
their lunar theories have first determined the mean longitude of the moon in 
functions of the true longitude as the independent variable, in the same manner 
as we have done in §4 of this chapter; and as the co-ordinates of the funda- 
mental differential equations employed by them have precisely the same signifi- 
cance as in this chapter, it would seem that the several results arrived at should 
be identically the. same, provided the methods of development employed were 
equally general and the calculations had been correctly made. 

La Place gives the two following equations for the determination of the 
mean longitude and the reciprocal of the radius vector, in functions of the 

8 
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true longitude as the independent variable, which are independent of the 
perturbations : 

nt = v — 2e(l — \f) sin {v — <o)+ f c^sin 2(v — w) — Je'sin B{v — a}) 

; (146) 



-=l + e2 + e(l + e^-i/)co8(v-ft>)+iey^cos(3v-a>~2Q) , ,^^^^ 

r \ ; (14b) 

+ -Jey^cos {v + (o — 2Q) 



while equations (89) and (105) will give, to terms of the same order, by neglecting 
the constant introduced by integration. 



nt = v — 2e sin (v — ft>) + f c^ sin 2(v — Q))—^e? sin 3 (v — fw) 

■ ; (147) 

+ \fQm2(v-Q)--iefsm{Sv-a}-2Si)-\f&m{v-Ba} + 2Q) " 



^=l-h^-{-e{l+e^coa{v-a))+iefcoQiSv-a)-2Q) . ,^^^^ 

r f . (148) 

- Jey2 co8(v + fli - 2Q) + ie>^ cos(v - 3a> + 2£J) 



Comparing now the values of nt given by equations (145) and (147), we notice 
that the coefficients of sin {v — w) in the two equations differ by the quantity Jey*, 
or by a term of the third order. We also see that La Place's equation contains 
a term of the third order depending on sin {v + w — 2Si), which does not appear 
in our equation ; while our equation contains a term of the third order depending 
on sin (v — 3a> H- 2Q), which does not appear in the equations of La Place or 
Plana except with a coefficient of the fifth order. These two last-mentioned 
differences arise directly from the introduction of the latitude of -the perigee into 
our analysis, an element that has never been explicitly considered in the analysis 
of the lunar and planetary theories, but which will be found to affect all terms of 
the third and higher orders in the planetary theories to the same extent as in the 
lunar theory. Indeed, the usual form of development is such that the argument 
of an inequality indicates the analytical order of magnitude of its coefficient, 
because the quantity (o is always connected with e, 2co with e^, etc. ; while S3 is 
always connected with y, 2Q> with y^, etc. ; so that the combination of the terms 
containing the sines and cosines of angles into which these quantities enter, pro- 
duce coefficients with powers and products of e and y in which the sum of the 
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indices is equal to the arithmetical sum of the coefficients of <d and S3 in the 
expression of the angle. But the introduction of the latitude ^o, of the perigee, 
brings w into the analysis unaccompanied by e, and in such manner that we have 
cos^o= {1 + >^sin^(a> — ^)}"i, which, being developed, produces the second 
mertiber of equation (77), an expression in which the order of the coefficient of 
any cosine is only one-half as great as the sum of the coefficients of w and Q> in 
the given angle. 

It is therefore easy to see that the product of two or more sines or cosines 
having coefficients whose order of magnitude has different relations to the con- 
stant terms. of the angles will completely vitiate the law in regard to the order of 
magnitude of any coefficient which would take place, provided the coefficients of 
all the sines or cosines bore the same relation to the constant terms of the angles. 
And hence the term depending on sin (v -f ai — 2SJ), of La Place's and Plana's 
analysis, disappears, and the term depending on sin (v — 3tt> + 2Q) comes in with 
a coefficient of the third order. 

If we now suppose that ft> = ^, in equations (145) and (147) they both reduce 
to the following : 

t? = n^ — (2e + J^)sin(v — ft>) + |^sin2(v — ft>) — (ie* + :J^ey*)sin3(i;--ft>) I nAG\ 
+ iy2sin2(v~Q) J 

In comparing equations (146) and (148), we notice that the coefficients 
of cos {v — w) differ by the term ^ef, while equation (146) has the term 
^ey^ cos {v + (0 — 2Q), instead of the two terms, — Jey^ cos {nt + <o — 2Si) 
-f ^ef cos {v — S(o+ 2Si)f in equation (148). However, if we suppose that a> = Q, 
they both reduce to the following : 

^ =l + ^ + e{l -he"- if) cos(v-a}) + iey^coa3(v- to). (150) 

r 

If we suppose, however, that v == tt>, in equations (146) and (148), they will 
become respectively 

- = 1 + e + c^ '\-e'''lef + lefcos2(o} - Q), (151) 

r 

and ^ = l + e + ^ + ^. (152) 



But the true elliptical value of - , when the moon is in perigee, is 

T 

= 1 + 6 + ^ + ^, which is the same as equation (152) ; and therefore 



r a(l-e2) 
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equation (151) is erroneous, unless oi = Si, in which case it becomes the same as 
equation (152), If we now take the differential coefficient of equation (146), we 
shall obtain 

a— =e(l+ e^-iy^ 8in(t; ~ai)+ fey* sin(3t; - a> ~2Q)+ J6j^sin(t;+ai-2a). (163) 

It is evident that this value of -- never becomes equal to nothing when v = ai, 

dv 

unless m is also equal to Q, or to Q + 90®, a conclusion already deduced from 

equation (16) of the Introduction. This defect, however, does not exist in equa- 

dir 
tion (146), for it gives — = 0, when v = ai, whatever be the relation between 

dv 
the quantities o) and Q. 



> 
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